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ABSTRACT 


Signal-processing theory for the TurboRogue receiver is presented. The signal form is 
traced from its formation at the GPS satellite, to the receiver antenna, and then through the 
various stages of the receiver, including extraction of phase and delay. The analysis treats the 
effects of ionosphere, troposphere, signal quantization, receiver components, and system noise, 
covering processing in both the “code mode” when the P code is not encrypted and in the “P- 
codeless mode” when the P code is encrypted. As a possible future improvement to the current 
analog front end, an example of a highly digital front end is analyzed. 



SYMBOL DEFINITIONS 


Roman 


A_,A 0 ,A + = 
A c = 

Am := 

a n 

A p (f) 

A x — 


early, prompt, and late correlation amplitudes, measured 

amplitude factor defined in Eq. (4.23) 

maximum correlation amplitude, Eq. (7.1) or Eq. (7.20) 

P-code chip value, ± 1 

P-code Fourier component, Eq. (2.6) 

amplitude factor defined by Eq. (5.7) 


C 


ci 

f 

f d 

fh 

fL 


= quantized counterrotation phasor, Eq. (4.2) 

= coefficient of fundamental in quantized counterrotation sinusoid 
= frequency 

= Doppler-shifted baseband frequency 
= nominal downconversion frequency 
= L-band carrier frequency 
= P-code chip rate = 1/T p 


G A (v d ) 

G(f) 

Gn(fd) 

G s (f) 

H (f,f') 

kB 

k p 

/ 

L c 

L s 


effective complex filter introduced by antenna 
composite bandpass amplitude, Eq. (3.12) 

effective aggregate complex filter introduced by receiver before sampler 
effective aggregate filter introduced by GPS satellite 

definition in Eq. (4.16) 

Boltzmann's constant 

constant in equation to convert from correlation amplitudes to residual delay 
lag offset relative to model delay 

SNR loss due to quantization of counter-rotation sums, Eq. (7.6) 

SNR loss due to sampling quantization, Eq. (7.7) 


v 



L x = model delay offsetting LI during cross-correlation of LI and L2 

N = (N s - l)/2 

N c = number of chips in accumulation 

N s = number of sample points in accumulation interval 

P(t) = P-code sequence 

Pn = average noise power entering the sampler, per quadrature component 

Pn = total noise power at RF at insertion 

Ps = average baseband signal power per quadrature component before receiver filtering 

qM = proportionality constant accounting for sampling quantization in low-SNR limit 

R cx = ratio of code and P-codeless SNRs, Eq. (7.26) 

R p (r) = time-averaged P-code autocorrelation function, Eq. (2.14) 

R p (t,r) = autocorrelation function for P code, Eq. (2.12) 

s = sample interval = (20.456)" 1 (Js 

S c = definition in Eq. (4.8) 

Si(t) = ideal unfiltered L-band signal, Eq. (2. 15) 

SNR 20 = "18-ms" voltage SNR 

SNR V = voltage SNR, Eq. (7.5) 

SNR X = maximum voltage SNR for LI *L2 correlation, Eq. (7.23) 

S p (f) = P-code power spectrum, Eq. (2. 1 1) 

S R (t) = signal received by the TR antenna, Eq. (2. 19) 

S T (t) = transmitted L-band signal, Eq. (2.18) 

S x = definition of Eq. (5.4) 

T = accumulation interval 

t = true time 

t = true time at interval center 
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T s 

U I (t R ) 


Vk 


V m (tR ,1) 

V(t) 

w 

x 

Xl,X 2 


Greek 

a f 

Acpi 

Wi 

R 

A(t) 

c?A 

CTn 

O u 

CT <t» 

X 

x' 


= (10.23) -1 |is = P-code chip interval 
= receiver time as indicated by sampler 
= system temperature 
= correlation sum, Eq. (4.5) 

= quantized baseband voltage representing signal after sampling, k = 1 for LI and 
k = 2 for L2 

= model signal generated by receiver, Eq. (4.1) 

= complex voltage representing signal at baseband, Eq. (3.2) 

= effective double-sided noise bandwidth 
= quantized value for quantity x 
= quantity x for LI and L2, respectively 


= loss in signal power due to receiver filtering 

= difference of instrumental phase between LI and L2, defined by Eq. (5.5) 
= bandpass average of Atpi 

= noise added to signal, at baseband before sampler 

= triangle function 
= 1 - cr amplitude error, Eq. (7 . 1 5) 

= RMS noise on either component of baseband signal before sampling 
= 1-G delay error, Eq. (7.19) or Eq. (7.20) 

= 1-G error on correlation sum 

= 1-G system-noise error in phase, Eq. (7.14) 

= combined delay, Eq. (2.20) 

= combined delay, Eq. (3.7) 

vii 



measured value for total delay, Eq. (6.9) 
aggregate group delay, Eq. (3.14) 
group delay due to ionosphere 
geometric delay from satellite to receiver 
composite instrumental group delay 
model group delay generated by receiver 

actual value for P2 group delay derived from P2-P1 group delay, Eq. (6.24) 

measured value for P2 group delay derived from P2-P1 group delay, Eq. (6.23) 

residual delay = difference of aggregate group delay and model delay, Eq. (4.24) 
and Eq. (5.9) 

error in receiver clock relative to true time, Eq. (3.5) 
measured value for residual delay, Eq. (6.7) 
error in GPS satellite clock relative to true time 
tropospheric delay 

actual value for phase delay, Eq. (6.6) 
measured value for phase delay, Eq. (6.5) 
aggregate phase delay, Eq. (3 . 1 6) 

actual value for L2 phase delay derived from L1*L2 phase, Eq. (6.19) 

measured value for L2 phase delay derived from L1*L2 phase 
measured value for total P2-P1 group delay, Eq. (6.20) 
actual value for total P2-P1 group delay, Eq. (6.21) 

Doppler-shifted frequency at RF, Eq. (3.3) 

actual value for total carrier phase, Eq. (6.4) 

measured value for total carrier phase, Eq. (6.3) or Eq. (6.1 1) 



9a 


= aggregate baseband carrier phase, Eq. (3 . 1 5) 

(p d = difference of aggregate baseband phase and model phase, Eq. (4. 10) or Eq. (5.8) 

(p e = ionospheric phase shift 

cp h = do wncon version phase, Eq. (3.1) 

(p hd = drift in downconversion phase relative to fht, Eq. (3.1) 

cpj = composite instrumental phase shift excluding tj term 

cpi = bandpass average of 9 i 

9L1 = measured value for total C/A carrier phase 

cpL 2 = actual value for total L2 carrier phase derived from L1*L2 carrier phase and C/A 

carrier phase, Eq. (6.17) 

9L2 = measured value for total L2 carrier phase derived from LI *L2 carrier phase and 

C/A carrier phase, Eq. (6.16) 

(p m = model carrier phase generated by receiver 

(p r = measured value for residual phase, Eq. (6.1) 

cp s = phase instability introduced by satellite, excluding frequency reference 

cp x = actual value for LI *L2 carrier phase, Eq. (6.14) 

<p x = measured value for LI *L2 carrier phase, Eq. (6. 1 3) 

y = "correction" term in 9 l 2 > Eq- (6.18) 

Other 

< > = ensemble average over both noise and code 

< >ri = ensemble average over noise, 77 

Op = ensemble average over possible code values 
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SECTION 1 


INTRODUCTION 

This report presents the nonrelativistic theory for signal processing in the TurboRogue 
(TR) receiver. The signal form is traced from the Global Positioning System (GPS) satellite to 
the receiver antenna, and then through the various stages in the receiver, including extraction of 
phase and delay. The analysis includes the effects of the ionosphere, the troposphere, discrete 
sampling, instrumentation, and system noise. The theory covers processing in the "code mode" 
when the P code is known and in the "P-codeless mode" when the P code is encrypted. 

Section 2 analyzes the GPS signal by performing a spectral decomposition and 
autocorrelation of the ideal P code and propagating the signal from the GPS spacecraft to the 
receiver, with the inclusion of both troposphere and ionosphere effects. Section 3 analyzes 
receiver front-end processing, tracing the signal from reception at RF to sampled form at 
baseband. Section 4 presents correlation theory for the code mode when the sampled signal is 
cross-correlated with a model signal generated by the receiver on the basis of feedback in phase 
and delay. Section 5 presents correlation theory for the P-codeless mode in which the LI and L2 
signals are cross-correlated. Section 6 describes extraction of the phase and delay observables in 
both the code mode and P-codeless mode and presents the corresponding theoretical expressions 
containing the various terms contributing to each observable. Section 7 analyzes the effects of 
system noise, including signal-to-noise ratios (SNRs), phase errors, and delay errors for both the 
code mode and the P-codeless mode. 

As an example of a possible future improvement to the current analog front end. 
Appendix A presents the theory for a particular design of a highly digital front end. Appendix B 
analyzes the effects of two-level and three-level quantization in the analog-to-digital conversion 
of the received signal and of three-level quantization of the counterrotation sinusoids. 

Appendix C presents noise correlation analysis, including correlation of noise between samples 
and correlation of noise on correlation sums with different lags. Appendix D analyzes the 
correlation of noise between the quadrature components of the baseband signal. Appendix E 
formulates the correlation function for the clear/acquisition (G/A) channel. 
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SECTION 2 


THE GPS SIGNAL 


2.1 Ideal P Codei 

The ideal, unfiltered P code, which is generated at a rate of f p = 10.23 MHz, consists of a 
pseudorandom sequence of chips, each with a value -1 or +1 and a width of T p = l/f p = 

(10.23)' 1 |is, as illustrated by plot (a) in Fig. 2-1. A particular P-code sequence can be expressed 
in its ideal form as 


P(t) = Xan B (t-t n ) (2.1) 


where the time argument t is true time, B(t) is the chip function shown in Fig. 2-2, and the 
coefficients a n are the chip values of -1 or +1. As indicated by Eq. (2.1), the chip centers fall at 
true-time values t n , which are equally spaced and given by 

t,=nT p -| (2.2) 

in which n is an integer whose range is sufficient to cover the time interval of interest. As 
implied in this time relation, the rising edge of a chip that falls at an integer-second boundary 
will be aligned with that boundary when the GPS clock has zero error with respect to true time. 
Satellite clock error will be accounted for below. 

The chip values a n will be approximated as random variables with zero mean: 

<a n>p = ° (2.3) 

and with no correlation between different chips: 

<aja n > p =8jn (2.4) 

where < > p indicates an ensemble average over possible chip values and 8j n is the Kronecker 
delta function. A given P-code sequence is generated for one week before the sequence is 
repeated. As developed below, the code sequences leaving the GPS satellites do not have the 
ideal form in Eq. (2.1) but are filtered to some degree. 

An expression similar to Eq. (2.1) also describes a C/A-code sequence with its chip 
width of (1.023)- 1 p s and a sequence repetition interval of a millisecond. With regard to 
correlation analysis, processing of the P and C/A data are very similar. As outlined in 
Appendix E, however, the short (1 -ms) repetition interval for the C/A sequences becomes a 
complication, but most aspects of the results remain the same. 


1 For an explanation of symbols used in this document, see page v, “Symbol Definitions.” 
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Figure 2-2. The P-Chip Function 


2.2 Spectral Analysis of the Ideal P Code 

A Fourier transform of a given P-code sequence, over the chips in an "accumulation 
interval" with center at true time t and length T, is given by 

J r t+T/2 

P(t) e‘ 27lift dt (2.5) 

t-T/2 

As assumed below, the accumulation interval will correspond to a correlation interval (e.g., 

18 ms) applied during data reduction. With the use of Eq. (2.1), the transform is given to 
excellent approximation by 


A p (f ) 


N c 

, sin(TtfTp) 
p TCfTp Zj 


a k e- 27lift “ 


( 2 . 6 ) 


where the sum over k applies only to chips within the transform interval and N c is the total 
number of chips in the interval. (A box will be used to enclose equations of particular 
significance.) This equation is not exact since the two end chips in the transform interval may 
not fall completely in the interval. This effect is negligible since the interval contains so many 
chips (e.g., about 180,000 in 18 ms). The frequency components will vary from correlation 
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interval to correlation interval as a function of the particular values for a^ (i.e., the particular 
P-code sequence within the interval), but will have zero mean over many intervals. The power 
spectrum, however, averages to a (sinx/x) 2 form, as shown below. 

Based on Eq. (2.6), the correlation between spectral components at different frequencies 
is given by 


<A;(f)A p (f’)> p = Tp 2 


sin(7tfT p ) sin( 7 tf T p ) 
rc fT P TCf'Tp 


Y <aj a k > p e 27li ( ft j- ft k) 
j.k 


By using the assumption in Eq. (2.4) and the definition of tk in Eq. (2.2), one obtains 


(2.7) 


<A£(f)Ap(f’)> p =T p 2 


sin(7tfT p ) sin(7tf T p ) 
^fT p jtf’Tp 


N c 

e*i(f'-f)T„Y 


k 


e 27ti(f-f')kT p 


( 2 . 8 ) 


One can show that the sum becomes 


e 7ti(f'-f)T p 'V e 2iti(f-f')kT p = e 27ii(f-f')t sin i N ctC(f-f )T p ] ^ g 

k sin[7c(f-f')T p ] 

where t is true time at the middle of the interval. This equation can be substituted in Eq. (2.8) to 
yield 


<Aj(f )A p (f , )> B = Tp e 2ici (f-f')i sm (rc f Tp) sm(7cf T p ) sin[N c n(f-f )T p ] ( 2 .10) 
P ^fTp 7tf'T p sin[7t(f -f )Tpj 

One can show (see Subsection 4,3) that the sinx/siny function can be accurately approximated by 
a sum of Dirac delta functions with a spacing of 10.23 MHz, provided N c is very large. Thus, the 
correlation between frequency components is substantial only when f - f 1 = n* 10.23 MHz, where 
n is an integer. 

The power spectrum is obtained by setting f = f , which gives 


S p (f ) = T p T 


sin(TcfTp) 


12 


7CfTp 


( 2 . 11 ) 


where the total accumulation time is given by T = N c T p Thus, on average, the power spectrum 
of the ideal P code is a (sinx/x) 2 function with a first zero crossing at (Tp)' 1 . 

2.3 Autocorrelation of the Ideal P Code 

The autocorrelation function of the ideal P code is defined by 
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( 2 . 12 ) 


R P (t,x) = <P(t)P(t+x)> p 

Autocorrelation of an example P-code sequence is illustrated in Fig. 2-1. Plot (b) offsets the 
P-code sequence in plot (a) by a selected delay x, while plot (c) shows the product of plot (a) and 
plot (b). Plot (d) shows the autocorrelation function obtained by averaging plot (c) over all 
possible P-code sequences. Plot (e) shows the particular "noise" sequence that separates plot (c) 
from plot (d). As implied, the "noise" sequence averages to zero over all possible P-code 
sequences. 

For 0 < x < T p , the autocorrelation function has a period T p as a function of time and a 
value of either 0 or 1 , as illustrated in plot (d) for a particular value of x. For x = 0, the 
autocorrelation function in plot (d) would be 1.0 at all times. For x > T p , the autocorrelation 
function would be zero. 

Averaging over a time interval containing a very large number of chips produces the 
time-averaged autocorrelation function 


R p (x) = 



R p (t,x) dt 


(2.13) 


Since the width of the positive pulses in plot (d) decreases as x increases according to Tp - x, the 
time-averaged autocorrelation function becomes 


R p (x)=A(x) 


(2.14) 


where A(x) is the triangle function shown in Fig. 2-3. The Fourier transform of a triangle 
function is (sinx/x) 2 function like that in Eq. (2.1 1), as one would expect. 

If all P-code frequency components were present in the signal (i.e., an ideal P code), this 
exact triangle function would describe the amplitude response produced by correlating the 
received GPS signal with an ideal model code. As shown below, filtering modifies the triangle 
response and smooths the sharp comers. 


2.4 The Received Signal 

This subsection develops the spectral form of the GPS signal arriving at the receiver, 
including geometric, tropospheric, and ionospheric effects. 

The form for the ideal unfiltered P code is given by Eq. (2.1). When modulated on an 
ideal carrier at L band, the ideal P signal becomes 

Si(t) = P(t)cos(27tf L t) (2.15) 

where fL is the L-band frequency. The L2 signal is constructed solely of a P signal of this form 
while the LI signal is a sum of such a P component and a similar C/A component. The C/A 
component is the same form as Eq. (2.15), but has a C/A pseudorandom sequence in place of P(t) 
and a carrier 90° out of phase with the P component (see Eq. (E.14) in Appendix E of this 
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document). The C/A component is transmitted with a power nominally 3 dB greater than the 
Ll-P component. 

The spectral decomposition of the P code can be written as 


P(t) = 



Ap(f)e 2jtift df 


(2.16) 


where the Fourier components are given in Eq. (2.6). This expression can be used to recast the 
ideal unfiltered L-band signal in Eq. (2.15) in the form 


S,(t) = A p (f)e 27ci(fL+f ) t df + c.c. (2.17) 


where c.c. denotes complex conjugate of the first term. The relation, A^(f ) = Ap(-f), which is 
valid for a real signal like P(t), has been used to obtain this formulation. In addition, an 
amplitude factor of 1/2 has been neglected since subsequent signal amplification at the spacecraft 
will lead to an overall gain factor as modeled below. 



Figure 2-3. The Time- Averaged Autocorrelation Function for the P Code 
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The signal leaving the spacecraft is filtered so that only the main lobe and a fraction of 
the first sidelobes remain. If G s (f) represents the effective aggregate (complex) filter and 
amplification introduced by the spacecraft, the transmitted signal as it leaves the spacecraft can 
be cast in the form 


S T (t)= I G s (f )A p (f )e 27ti [( f L+f)(t+ts)+9s] df + c.c. (2.18) 


where cp s represents phase instability in the satellite frequency system, excluding the frequency 
reference. The spacecraft filter is reportedly a six-pole Chebychev with 0.01-dB ripple and a 
3-dB double-sided bandwidth of about 36 MHz. A clock offset x s has been included to represent 
the error in the satellite clock relative to true time. (Satellite time = true time + t s .) For 
simplicity, this expression neglects the multiplicative data-bit factor contained in the actual GPS 
signal. The following analysis assumes that a correlation interval is less than a data-bit width 
and that data-bit synchronization has been established. In this case, the data-bit factor would 
only change the analysis by a multiplicative sign factor throughout the equations. 

In transit to the ground, the signal is retarded by a geometric delay, t g , and a tropospheric 
delay, and is advanced in phase by an ionospheric phase shift, cp e . After incorporation of these 
effects, the received signal arriving at the receiver antenna at true time t can be written in the 
form 


S R (t) 



Gs(f)A p (f ) e 27ti[(fL+f)(t-T) + ( p s + ( p e ] df + c.c. 


(2.19) 


where the combined delay is given by 


T=T g +X t -X s (2.20) 

Decreases in signal strength in transit are absorbed in G s . The RF signal in Eq. (2.19) is 
collected by the receiver antenna for processing. 
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SECTION 3 


FRONT-END PROCESSING 


This section analyzes the steps in the TurboRogue front end that take the signal from 
reception at RF down to its form after baseband sampling. The current analog front end is 
analyzed in this section and an example of a highly digital front end is analyzed in Appendix A 
of this document. 


3.1 Functional Description of the Analog Front End 

The receiver antenna collects the GPS signal as input to the analog processing, as shown 
in Fig. 3-1. The output from the antenna is passed through an RF bandpass filter spanning 1075 
to 1725 MHz to eliminate out-of-band RF interference (RFI). The signal is then amplified by a 
gallium-arsenide field-effect transistor (GaAs FET). System temperature, which includes all 
noise sources, is nominally 200 K for LI and 225 K for L2. The resulting signal is power 
divided into LI and L2 branches and each branch is downconverted to baseband with a 
quadrature mixer in a double-sideband mode with a fixed-frequency local-oscillator (LO) signal. 
The four resulting baseband signals are each passed through a low-pass seven-pole Butterworth 
filter with a single-sided bandwidth of 9.6 MHz, as plotted in Fig. 3-2. The four filtered signals 
are then each hard limited and "sign" sampled (1-bit quantization) at a rate of 20.456 MS/s. 

(MS/s will be used to denote 10 6 samples per second.) The same LI samples are processed by 
the digital signal processor to obtain both the C/A and the P observables. 

To carry out this processing, three frequencies are required: the "sample-clock" rate at 
f s = 20.456 MHz and two LO (downconversion) signals, fhi = 1575.1 12 MHz for LI and fh 2 = 
1227.36 MHz for L2. The frequency-and-timing subsystem generates 20.456 MHz from a 
5-MHz reference and then generates the two LOs by simple multiplication of the sample-clock 
rate. That is, the LI LO is generated by the multiplication 77f s (which equals fy - 308 kHz) and 
L2 by 60f s (which equals fy - 240 kHz, where f u and f L2 are the zero-Doppler carriers at RF). 

Care must be taken in selecting the sampling frequencies and downconversion 
frequencies. The sampling frequency must be essentially incommensurate with the code chip 
rate (i.e., 20.456 MHz vs. 10.23 MHz) so that discrete-sampling errors will be negligible. 
Theoretical verification of the adequacy of 20.456 MHz is presented below in Subsection 4.3. 

To produce positive baseband carrier frequencies, each downconversion frequency is offset from 
its respective carrier (by 308 kHz for LI and 240 kHz for L2). The large resulting positive 
baseband carrier frequencies reduce to negligible levels the errors caused by the three-level 
quantization of counterrotation sinusoids applied in subsequent processing. Subappendix B.2 
estimates the errors introduced by carrier counterrotation. 


3.2 Mathematical Model for the Analog Front End 

As described in the preceding subsection, the receiver front end passes the received signal 
represented by Eq. 2.19 through amplifiers and filters and downconverts from L band to 
baseband. This subsection models the composite effect of the antenna and front-end steps as a 
linear process that can be formulated as a single overall complex system filter and a single 
effective phase and frequency shift representing the downconversion process. 
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LI: COS 
20.456 MS/s 


LI: SIN 
20.456 MS/s 


L2: COS 
20.456 MS/s 


L2: SIN 
20.456 MS/s 


Figure 3-1. Block Diagram for the Analog Front End 
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Do wncon version can be modeled in the following way: Let the do wncon version phase at 
the mixer be denoted by <ph. Because the frequency reference of the receiver is very stable, 
downconversion phase can be represented by a time function with the ostensibly linear form 

<Ph = fht + <Phd (3.1) 

where fh is the nominal downconversion frequency and <Phc is a slowly varying drift relative to 
the nominal linear function associated with fj,. The quadrature mixer multiplies the (real) 
received RF voltage by cos cph and by sin %, a process that can be modeled in complex form as a 
multiplication by e- 2lli< Ph. When the expression in Eq. (2.19) is multiplied by e- 2ui( Ph ; the 
multiplication and integral can be interchanged in both the first term and the c.c. term so that the 
multiplication is applied separately to each frequency component. For each frequency 
component, one obtains a small difference frequency from the first term and a large sum 
frequency from the c.c. term. Since the large sum-frequency term is eliminated by filters 
following the mixer, only the difference frequency term must be retained in the analysis. 

The preceding e- 27ti< P>> model for downconversion from RF assumes perfect 90° separation 
of the sine and cosine mixing signals. In practice, the analog front end will produce sinusoids 
that deviate from perfect 90° separation. This offset from 90° will be called the quadrature- 
separation error. A quadrature-separation error introduces incoherency between the two final 
complex correlation sums obtained by processing the two quadrature components. Relative to a 
perfect amplitude of 1.0, the lowered output amplitude would be cos (A0/2), where A0 is the 
quadrature-separation error in radians. When A0 is 16°, the voltage signal-to-noise ratio (SNR) 
is reduced by about 1%. According to D. J. Spitzmesser (Spitzmesser, 1994), the quadrature- 
separation error for the TurboRogue receiver is about 5° or less, which would cause an amplitude 
loss of about 0.1% or less. This result indicates that the quadrature-separation error for TR is of 
little ultimate importance and it will therefore be ignored in the following analysis. 

Antenna gain and possible antenna filtering, phase shifting and delays will be modeled by 
a complex antenna filter, G A , that is a function of RF frequency. Similarly, all amplification, 
filtering, phase shifts, and delays due to the front end will be modeled with a single complex 
filter, G r . This front-end filter consists of a product of the filters at successive stages as well as 
phasors representing other phase shifts and delays. Frequency-component identification for the 
front-end filter, however, will be specified by means of a Doppler-shifted baseband frequency, 
and, as a consequence, a translation is required to make an association with the corresponding RF 
frequency preceding downconversion. This convention is motivated by the fact that the 
dominant filter is the "narrow" Butterworth filter applied at baseband. 

As indicated in Fig. 3-1, each L-band channel has two quadrature components and each 
of these components has a separate Butterworth filter. According to Spitzmesser (Spitzmesser, 
1994), the four Butterworth filters are not identical but are matched to within about 7% with 
respect to bandwidth. To account for the differences, a detailed analysis would provide a 
separate function for each of these filters. One can show, however, that the final effective filter 
for a given L-band channel is a straight average of the two Butterworth filters for that channel. 
Thus, for simplicity in analysis without sacrifice of accuracy, it will be assumed at this point that 
the two Butterworth filters for a given L-band channel are identical, each with a response equal 
to the average of the two actual filters. 

Another potential concern is how well delay and phase shifts are matched between sine 
and cosine components in their paths from the RF branching point to baseband sampling. 
According to Spitzmesser (Spitzmesser, 1994), the path length for each component is only about 
20 cm, and the two paths are matched to about 3 cm (0. 1 ns). The relative path difference within 
the quadrature mixer before mixing is about 0. 1 cm or less, and the corresponding relative phase 
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shift is about 0.1 cm/20 cm = 0.005 cycles (1.8°) The relative phase shift after mixing is about 
300 kHz * 0. 1 ns = 30 jlcycles. These relative phase and relative delay values are too small to 
cause significant loss of amplitude. Observed delay will contain the average of the two 
component delays and this implicit average delay, combined with other delays, will be either 
calibrated, differenced out, or estimated in subsequent processing. Thus, none of these effects is 
of ultimate importance and for this reason, none will be separately modeled. 

The preceding discussion indicates that the antenna and front end can be modeled by 
inserting the antenna filter and the composite receiver filter into the integrand of Eq. (2.19) and 
multiplying by e‘ 27ti( P h . Since filtering in the front end removes the sum-frequency term 
resulting from the c.c. term, the complex voltage (real and imaginary) at baseband becomes 


V(t)= G R (f d )GA(v d )Gs(f)A p (f)e 27ci [(fL+f)(t-T)+cp s +(Pe-9h]df+Ti (3.2) 


where the Doppler-shifted frequency at RF corresponding to f, the baseband variable of 
integration, is given by 


v d = (f L +f)(l-x)+ <p s + cpe ( 3 - 3 ) 

in which the top dot denotes time derivative. The corresponding baseband frequency is given by 

fd= Vd- fh (3.4) 

The complex term T) represents all noise, after the downconversion and filtering operations, 
including the noise picked up by the antenna and the noise added by the receiver. Almost all of 
the noise contributed by the receiver is added at RF. Analyses of the effects of the noise term are 
presented in Section 7, Appendix C, and Appendix D of this document. 


3.3 Mathematical Model for the Sampled Signal 

The baseband signal theoretically represented by Eq. (3.2) is hard limited and sign 
sampled at a sample rate of 20.456 MHz. Timing for the sampler is derived from the receiver 
clock, which will have an offset from true time. True time t and receiver time t R (as indicated by 
the sampler) will be related by 


t R — t + Tr (3.5) 

where x R is the error in the receiver clock. A quantized complex sample value at time t R will be 
denoted by V(t R ). For a particular time t R , the theoretical form for the voltage Y is obtained by 
substituting t R - t R in place of true time t in the right-hand side of Eq. (3.2). 

As discussed in Appendix B for low SNR signals, the sampled voltage after averaging 
over noise is approximately proportional to the signal component divided by the root-mean- 
square (RMS) noise. Thus, when the complex voltage in Eq. (3.2) is subjected to two-level 
quantization in each quadrature component and then to an ensemble average over noise, rj, the 
theoretical form for the averaged complex voltage becomes 
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<V(t R )> 11 = ^ 

0^ 



G R (f d )G A (v c i)Gs(f )Ap(f )e 2jti KfL+f)(tR-x')+cp s + (p e - <p,0 df 


(3.6) 


where the tilde on V denotes amplitude quantization, < ^ denotes average over noise and On is 
RMS noise on either quadrature component. (It will be assumed that RMS noise is the same for 
the two components. The fact that this assumption is only approximately valid does not change 
the results of the following analysis, since only single-sample SNR matters in 1-bit sampling, 
and single-sample SNR is the same for the two components.) Delay has been redefined to 
incorporate the sampler clock error that enters when true time is replaced with sampler time: 

T' = T + t R = tg + tt - t s + Tr (37 

As derived in Appendix B, the proportionality factor qM accounts for sampling quantization and 
is given by 



(3.8) 


for two-level quantization. 


3.4 Reformulation of the Sampled Signal 

Subsequent analysis can be clarified by explicitly defining and collecting terms of like 
nature in the signal component. First expand the ionosphere phase about the carrier frequency 

Cpe(fL+f) ~ (Pe(fL) + -Jr f ( 3 -9) 

which becomes 

(fe(f L +f) = (pe(f L )-'Cef (3.10) 

where t e is the group delay due to the ionosphere. Next, reformulate and consolidate the 
bandpass filters 


G R (fd) G A (Vd) G s (f ) = G(f)e-2*i(f*.-<P.) 


(3.11) 


where the composite bandpass amplitude is defined as 

G(f) = |G R (f d )G A (v d )G s (f)| (3.12] 

and Ti is the overall instrumental group delay, including spacecraft, antenna, and receiver. The 
frequency-dependent term (pi is the instrumental phase effect that remains after removing the 
linear phase-versus-frequency function associated with ij. Note that the bandpass amplitude, the 
instrumental phase, and the effective instrumental delay can potentially depend on Doppler shift, 
as indicated in Eq. (3.1 1) by the Doppler-shifted arguments, fd and Vj. 
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The effective delay introduced by the Butterworth fdter described in Subsection 3.1 is 
about 77 ns for the P code. This delay is the approximate theoretical value for the Butterworth 
filter itself and does not include other instrumental delays introduced by the front end. 

When Eqs. (3.10) and (3.1 1) are substituted in Eq. (3.6), the noise-averaged baseband 
signal becomes 


<V(t R )> 11 =^e 2lti( Pa G(f)A p (f)e 2lti [f(tR-T a )+(p,] ( if (3.13) 

a n I 

* -oo 

where aggregate group delay is defined by 


X a a X g + X t - T s + Tr + X e + Xi 


(3.14) 


and aggregate baseband carrier phase by 


tPa =fc(t R -T(p) - (Ph 


(3.15) 


in which aggregate phase delay is given by 

X^ s Xg + T t - X s +X R - f£.'[<ps + (pe(fL)] ( - 3 ‘ 16 ' ) 

All time-varying quantities are evaluated at sample time, t R . 

Equation (3. 13) indicates that the baseband signal consists of a product of a carrier factor 
with a rate of fi/l - xip) - fh and a code factor involving an integral over the Fourier components 
of the P code. This formulation has representative variables for the major effects encountered in 
processing GPS signals. Note, however, that a bandpass-averaged instrumental phase shift 
emerges from the integral, and, furthermore, the best definition of an effective instrumental delay 
can be application dependent. 
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SECTION 4 


CORRELATION IN THE CODE MODE 


When the P code is known, the receiver correlates the baseband samples with a model 
signal consisting of the appropriate P-code sequence, offset by a model delay, and a carrier 
phasor incorporating a model phase. Model phase and model delay are supplied in real time as 
feedback by delay- and phase-locked loops operating on correlation output. This report will not 
treat feedback theory, but will assume that model phase and model delay are adequately provided 
by loops operating on correlation output. Correlation analysis for the C/A code, which is slightly 
different than P-code analysis, is presented in Appendix E. 


4.1 Model Signal 

The model signal generated by the receiver at receiver (sampler) time t R has the form 


Vm(tR>0 — P(tR-T m -Is) C((p m ) (4.1) 

where P is the ideal P-code function in Eq. (2.1), x m is the feedback model for total group delay 
at time t R , l is a lag offset relative to the model delay, s is the sample interval of (20.456)' 1 (is, 
and cp m is the feedback model for phase at time tR. Note that the receiver generates the code 
sequence as a function of t R , as though t R were true time. The model signal is generated for each 
of three lags (-1, 0, +1) in the TurboRogue receiver. As a function of time, the generated model 
code consists of the sequence of -Is and -t-ls specific to the satellite being tracked, but with 
effectively instantaneous transitions between chips (i.e., with no filtering), as implied by 
Eqs. (2.1) and (4.1). Further, the P-code sequence is generated with a delay that is effectively 
exactly equal to the model delay generated by the feedback loops. Over a correlation interval of 
18 ms, the receiver generates both the model delay and the model phase as linear time functions. 

The model signal contains both sine and cosine components in order to provide _ 
quadrature counterrotation of the carrier of the sampled signal. The complex quantity, C , is a 
quantized counterrotation phasor defined by 

C(<p m ) = cos(<p m ) - i sin((p m ) (4.2) 

where the tildes indicate three-level quantization, as explained in Appendix B. The quantized 
phasor can be expanded in terms of the harmonics of the fundamental in the form 

C(tpm) = cie- 2ni( Pm + harmonics (4.3) 

where ci = 1.176 is the coefficient of the fundamental for the chosen three-level quantization. 

Remembering that P(t) is real, one can take the complex conjugate of Eq. (2.16), 
substitute the result in Eq. (4.1), and obtain a useful form for the receiver model in terms of 
frequency components: 
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(4.4) 


V m (t R ,/)=C((Pm) A;(f')e-2^f , (tR-^-/s) d f 


where Ap(f) represents the Fourier components in Eq. (2.6) pertaining to the particular code 
sequence within the correlation interval. Note that the receiver time argument in Eq. (4.1) has 
been substituted directly in Eq. (2.16) without including the difference, x R , between receiver time 
and true time. This procedure is based on the fact that the receiver generates its P code exactly as 
represented in Eqs. (2.1) and (2.16) but on the basis of its own time scale: receiver time. No 
bandpass filtering is included in this equation because the receiver generates a P-code sequence 
that is essentially an ideal sequence. 


4.2 Correlation Sum 

At each sample point, the receiver generates the model signal in Eq. (4.1) and multiplies 
it by the sampled signal. For each lag value, the resulting product is accumulated for all points 
within a correlation interval to obtain a correlation sum, computed as 

N 

u/(t R ) = X V(t R +ks) V m (t R +ks, /) (4.5) 

k = -N 

where V is the complex sampled signal, t R is now the receiver time at the center of the sum 
interval, and the index k covers the sample points within the correlation interval. The correlation 
interval, whose length (integration time) will be denoted by T, is 1 8-ms long and contains about 
360,000 sample points for the current implementation of the TurboRogue receiver. In Eq. (4.5), 
it is assumed without significant loss of generality that the interval contains an odd number of 
samples with the total number given by N s = 2N + 1. The integration time is related to total 
number of samples by T = N s s. Since the sampled signal and model signal are complex, the 
correlation sums are complex. 


4.3 Correlation Function 

The correlation function is obtained by taking an ensemble average of the correlation 
sums in Eq. (4.5) with respect to all possible noise patterns and P-code sequences. An ensemble 
average over noise gives 


N 

<u;(t R )> T1 = X <M[t R +ks)> T1 V m (t R +ks, /) (4.6) 

k = -N 

where < >t, denotes the noise average. On the right side, the noise average has been 
interchanged with the sum and applied to the only factor containing noise, the sampled signal. 
After substitution of Eqs. (3.13), (4.3), and (4.4), and application of an ensemble average over 
code, Eq. (4.6) becomes 
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CO 


( 4 . 7 ) 


<U/(t R )>= Cl 


G(f) <Ap(f) A^(f')> p S c df df ' 


where 

N 

S c =^ e 27ti[(pa-<pm + (f-f')(tR+ks)-fT a +f , (tm+/s) + <pi] (4.8) 

k=-N 

in which the phase and delay functions are evaluated at time t R + ks. The unlabeled brackets < > 
denote ensemble average over both noise and code, while < > p denotes ensemble average only 
over code. To reach this form, the sum and integrals have been interchanged and only the 
Fourier components within the integral have been treated as random variables subject to the 
ensemble average over code. Further, only the fundamental of the counterrotation phasor in 
Eq. (4.3) has been retained under the assumption that the higher harmonics are filtered out by the 
sum over the correlation interval. 

Counterrotated phase in Eq. (4.8) can be expanded about the center time to yield 

(pa(t R +ks) - qj^tn+ks) ~ (pd+ 9 dks (4.9) 

where the phase difference at interval center t R is defined by 


<Pd= 9a(t R ) - <Pm(t R ) (4.10) 

and where <Pd is the feedback error for phase rate across the interval (i.e., the last equation but 
with top dots). Similarly, the delay terms in Eq. (4.8) can be expanded to yield 

fx a (t R +ks) - f 'x m (t R +ks)«f [t a (t R )+x a ks]-f ' [T m (t R )+i; m ks] (4.11) 

which becomes 

fx a (t R +ks) -f 'x m (t R +ks) «f-r a (t R ) - f^t*) + (f-f')x a ks (4.12) 

where, in the last equation, it has been assumed that the model delay rate very accurately tracks 
actual delay rate. One can show that this approximation is acceptable for even the worst-case 
errors in delay rate (<5 ns/s). 

When Eqs. (4.9) and (4.12) are substituted in Eq. (4.8), one obtains 

N 

g c ~ e 27li[<pd+(f-f')t R -f't a +f'(tm+/s) + (pi] ^ e 27ti[(f-f')(l-Ta)+<Pd]ks (4.13) 

k=-N 

where T a and t m are now evaluated at interval center. The sum in this equation can be evaluated 
in a manner similar to the derivation of Eq. (2.9), which yields 
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S c = e 27ti [ < Pd + ( f - f ’) tR- a+ f '(t m + 1 s) + 9,] sin|N s n[(f-f ) .U~ V +< ftd s ] 

sin[7t[(f-f')( l-t a )+q) ( j]s] 


(4.14) 


where N s is the total number of sample points in the interval. 

When Eq.(4.14) and the expectation value in Eq. (2.10) are substituted in Eq. (4.7), one obtains 

foo f oo 

<uKt R )> = ci^T^e 2lti< Pd I I H(f,f')e 27ti C(f-f'HR-fta+f'(xm+/s)+<pi]dfdf' (4.15) 

J - oo c 


where the two-dimensional "filter" function is defined as 

H(f f'^_ G (f) sin (tt fT p) sin(ttf'Tp) sin[N c 7t(f-f')T p ] sin[N s 7t[(f-f')(l-T a )-i-(i) d ]s] ^ ^ 

wf’Tp sin[jt(f-f')T p ] sin[jt[(f-f’)(l-'i^)+9 ( j]s] 

Equation (3.5) has been used to replace t R - 1 with t R . 

In order to schematically illustrate the nature of the "sinx/siny" functions in Eq. (4.16), 
Fig. 4-1 plots Eq. 4.16 as a function of f - f for an unrealistically small correlation interval 
(N s = 22) and for an excessive sample rate of 22.456 MS/s. These plots show that both of the 
two sinx/siny functions are "comb" functions consisting of a sequence of equally spaced peaks, 
with a peak spacing of 10.23 MHz for the first sinx/siny function and about 22.456 MHz for the 
second. Peak height is given by the number of chips (N c = 10) for the first sinx/siny function and 
by the number of samples (N s = 22) for the second. For both, the width is T/T, where T is the 
length of the correlation interval. For an actual correlation interval (e.g., T = 18 ms and a P 
sample rate of 20.456 MHz), these peaks become extremely high (e.g., about 10 5 larger than 
intervening peaks) and extremely narrow relative to the chip rate and sample rate (e.g, 55 Hz 
versus 10.23 MHz). For the example shown in Fig. 4-1, note that only the "origin" peaks at 
f-f ’= 0 line up between the two sinx/siny functions in Eq. (4.16), while the other peaks do not 
line up. This absence of peak coincidence away from f - f 1 = 0 , which is a consequence of the 
effective incommensurability between 10.23 Mhz and 22.456 MS/s, leaves only the origin peaks 
to make a significant contribution to the integral. If only the origin peaks are significant, the 
system has been successfully designed so that discrete sampling effects are negligible. 

For the actual sample rate of 20.456 MS/s, one can also show that only the origin peaks 
are significant. For zero Doppler, the commensurability ratio for the sample rate and chip rate 
reduces to 


f p_ s _ 10230000 _ 5115 (A 1T . 

f s T p 20456000 10228 

after common factors have been removed. Given this result, one can easily show for zero 
Doppler that the first coincidence of peaks for If -f I above zero occurs at 

|f — f ’| = 5115* 20.456MHz = 10228 * 10.23MHz (4.18) 


or equivalently 


|f- f ' | = 104.63244 GHz 


(4.19) 
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For | f - f ' | to be this large, the magnitude of either f or f ' must be greater than about 52 GHz. 
Because of the filtering effect of G and the two sinx/x functions at high frequencies in Eq. (4.16), 
this particular coincidence peak and other coincidence peaks at even higher multiples are reduced 
to negligible influence in the integral. For example, for a first coincidence peak with f ' = 

52 GHz, the sinx/x filter supplies an attenuation of about 0.00006 (-84 dB). One can show that 
an amplitude error of 0.00006 leads to a delay error on the order of 1 mm. 

Inserting nonzero values of Doppler (represented by delay rate, x a ) will improve effective 
incommensurability. Delay rate for ground-based receivers, for example, is less than 3 Jis/s. 

One can show that a Doppler effect this size will shift the peaks of the sinx/siny function by an 
amount too small to align any peaks that were not already aligned for zero Doppler. More 
importantly, Doppler shifting will dealign peaks that would have been coincident (but are not at 
the origin) for zero Doppler and thereby reduce commensurability effects to even smaller levels. 
On the basis of Eq.(4.16), one can easily show that Doppler reduces the amplitude of a given 
product of coincident peaks by a factor of about 7t Af x a Ti or more, where Ti is the integration 
interval and where Af = |f - f '| is the frequency difference at the peak. For the very small 
Doppler value of 1.5 Hz (x a = 1 ns/s) and the short integration time of 1 s, for example, the 
reduction factor is equal to 327 at the first coincident peak at 104 GHz. Thus, even for very 
small Doppler values, the commensurability error is greatly reduced by the Doppler effect. 

The preceding analysis indicates that commensurability of sample rate and chip rate will 
introduce delay errors on the order of 1 mm for the worst case of zero Doppler and even smaller 
delay errors when Doppler is included. 

When in lock, the rate-error term <p d in Eq. (4. 16) is not large enough to cause peaks to 
completely dealign, but can cause amplitude loss. This loss is accounted for below in the final 
form for the correlation function. 

When the number of sample points N s is very large, and when the sample rate and the 
chip rate are essentially incommensurate, as indicated, the above considerations show that the 
product of the two sinx/siny functions reduces to a spike function at f-f ' = 0. By approximating 
the two sinx/siny functions by sinx/x functions near f-f =0 and approximating the first of these 
two sinx/x functions by a Dirac delta function (assuming N c is very large), one can show that this 
product can be accurately approximated, within the significant ranges of the two sinx/x filters in 
Eq. (4.16), by 

sin[N c rc(f-f' )T p ] sin[N s 7t[(f-f )(1 -x a )-Kpd]s] ^ sin(N s 7t(p d s) g 
sin[7t(f-f )T p ] sin^Kf-f 1 )(1 -x a )+(pd]s] t p N s 7ttp d s 


where 5(f - f') is the Dirac delta function, given in this case by the limit 


8(f-f')= TpLIM 


sin[N c 7t (f-f 1 ) T p ] 
7t(f - f ') Tp 


N c ->o° 


(4.21) 


When Eqs. (4.16) and (4.20) are substituted in Eq. (4.15), the correlation function 
becomes 
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(4.22) 


<u/(t R )>= A c e 2lti( P<i 


r 

1 nff 

sin(rcfTp) 

I ; 

J-oo 

7tfTp 


e -2ici[f(Xr-/s)-<pi] df 


where the amplitude factor is defined by 


A c — N s TpCj 


q M sin(7t(pdT) 

TccpdT 


(4.23) 


in which the expression T = N s s has been used to simplify on the basis of correlation sum length. 
In this expression, difference phase, <Pd, is defined by Eq. (4.10) and residual delay by 


T r = T a (tlO - 1m(tR) (4.24) 

with both time tagged at the center of the correlation interval. Note that an additional phase shift 
can emerge from the integral in Eq. (4.22) to contribute to phase. This result for the cross- 
correlation function is very similar to the form that would be obtained in analog processing and 
therefore indicates successful removal of discrete-sampling effects. 


4.4 Discussion of the Correlation Function 

As one would expect, correlation amplitude Eq. (4.22) increases in proportion to N s , the 
number of sample points in the correlation sum. As discussed in Appendix B of this document, 
the multiplicative factor q M and the RMS noise On account for sampling quantization. The 
factor ci represents the amplitude increase due to tnree-level quantization of the counterrotation 
sinusoids. A sinx/x multiplicative factor models the amplitude loss due to feedback error in 
carrier-phase rate across the correlation interval. Rate feedback is accurate enough to keep 
amplitude loss resulting from this effect at an acceptable level. (Typically, this loss is negligible, 
but can reach 0.2 dB under worst-case dynamics.) A multiplicative phasor based on (Pd accounts 
for most of the feedback error in carrier phase at interval center. (A small additional carrier- 
phase effect can emerge from the correlation integral to contribute to the model for residual 
phase.) Typically, tracking error, including the contribution from the integral, is very small 
(< 0.02 cycle) but might reach 0.1 cycle under extreme dynamics. 

The integral across the passband yields a complex correlation response that is a function 
of model delay x m and lag l and provides a measure of the alignment of the model P code with 
actual P code. The integrand contains the system filter in Eq. (3.12) resulting from spacecraft, 
antenna and receiver, and a (sinx/x) 2 function representing the average power spectrum of the 
pure P code. This integral peaks in magnitude when the model delay, including lag offset, 
closely approximates the actual delay. After the correlation sums have been obtained from the 
correlation process, residual delay, can be estimated through analysis of the variation of 
amplitude with lag, as discussed in the next section. Typically, feedback is accurate enough to 
keep the model P code aligned with the incoming P code to better than 0.001 P-chip so that a 
given lag (the "prompt") produces an amplitude very near the peak of the correlation function. 
Correlation amplitude for Ll-P is plotted in Fig. 4-2(a) as a function of model delay on the basis 
of the bandpass function for the current analog front end described in Subsection 3.1. The model 
delay in these plots is referenced to true delay. (That is, a model delay of zero corresponds 
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Simulated with (sinx/x) 



Figure 4-2. P-Code Correlation Function for the TurboRogue Receiver 
with the Analog Front End 






approximately to zero residual delay.) Asymmetry in the amplitude-versus-delay response in this 
example is caused by the nonlinear phase- versus-frequency response of the receiver's baseband 
Butterworth filter. If the system filter passed all the frequency components of an ideal P code, 
the amplitude-versus-delay response would have the exact triangle shape shown in Fig. 2-3. 
Filtering rounds the sharp comers and broadens the response to the form shown in Fig. 4-2(a). 

Fig. 4-2(b) presents the corresponding phase-versus-delay response, which is nearly constant 
around zero residual delay. For purposes of illustration, the phase at zero delay in this plot has 
been offset from zero. In practice, this phase offset does not appear in the measured correlation 
sums corresponding to Eq. (4.22) since the phase-locked loop removes it along with all other 
phase offsets. Fig. 4-2(c), which presents the derivative of P correlation amplitude with respect 
to delay, is used in Subsection 6.1.2. 

For comparison. Fig. 4-3(a) plots C/A correlation amplitude as a function of model delay, 
given a pure (sinx/x) 2 power spectrum for the C/A code. (As discussed in Appendix E, code 
self-noise can cause significant deviations from this (sinx/x) 2 spectrum.) The C/A response is 
almost a perfect triangle, since the Butterworth filter is so broad relative to the C/A spectrum. 
Slight rounding near the peak and slight ringing near +1 c-chip is barely discernible. Fig. 4-3(b) 
presents the corresponding phase-versus-delay response. Similar to the P result, the phase at zero 
residual delay has been offset from zero. Note the slight ringing caused by the loss of sidelobes 
above the tenth sidelobe as a result of the Butterworth filter. Fig. 4-3(c), which presents the 
derivative of C/A correlation amplitude with respect to delay, is analogous to Fig. 4-2(c). 
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SECTION 5 


CROSS-CORRELATION IN THE P-CODELESS MODE 

When the P code is encrypted, the Ll-P and L2-P channels can be processed in a "P- 
codeless mode" that cross-correlates the two channels and thereby takes advantage of the fact the 
P-code sequence is the same on the two channels. In the following analysis, P-codeless 
operations will be denoted by L1*L2. For each sample point, the sampled voltages from the two 
channels are multiplied together, with one channel offset by a feedback model delay. The carrier 
of the product is counterrotated on the basis of feedback model phase, and the result is summed 
over a correlation interval of 18 ms. The resulting correlation sums are analyzed to generate the 
next values for phase and delay feedback and to extract the difference phase and difference delay 
between the Ll-P and L2-P signals. The primary disadvantage of operating in the P-codeless 
mode is that SNR is lower and output observables are noisier. 


5.1 Correlation Sum 

For each lag value, the receiver generates the following cross-correlation sum over a 
correlation interval of length T: 


N _ 

u/(t R )= X Vi(t R + ks-LxS-/s) vJ(t R +ks)C(cp m ) (5.1) 

k=-N 

where V] and V 2 are the (complex) sampled voltages for the LI and L2 channels, respectively; C 
is the counterrotation phasor in Eq. (4.2); <p m is model phase at time t R + ks; L x is the model 
delay fed back by the tracking loop; and / is one of the three correlation lags (-1,0, +1). Again, 
the total number of samples in the sum is odd and given by N s = 2N + 1 ; the correlation time is 
given by T = N s s; the time tag t R is at the center of the sum interval; and the index k covers the 
points in the interval. A theoretical expression for the noise- averaged value for the LI or L2 
baseband signal can be obtained from Eq. (3.13) by attaching the appropriate channel subscript to 
each channel-specific symbol. 

The delay difference between the LI and L2 channels is primarily due to the ionosphere 
and is therefore a relatively small, slowly varying difference. Based on analysis of the 
amplitude-versus-lag dependence of cross-correlation sums, the receiver feeds back a slowly 
varying model delay, L x , that is needed to place the maximum correlation response at the prompt 
lag. Because the L2 channel experiences a larger ionospheric delay than the LI channel, the 
model delay is applied to the LI channel so that lagging will be in the right direction to bring the 
two channels into delay alignment. Unlike the model delay applied to the code in the code mode, 
the model delay applied to LI in the P-codeless mode can only be changed in steps of a sample 
interval, which is given by s = (20.456 MHz)' 1 . Since the ionospheric delay difference between 
LI and L 2 is typically 100 ns or less, L x is typically 0 , 1 , or 2 (in units of sample interval, s). 

As shown below, phase of the V i*V 2 product is equal to LI phase minus L2 phase. The 
phase-rate difference between LI and L2 is equal to an LO term of 68 kHz plus a Doppler term 
that is about 1 kHz or less in magnitude for earth-fixed receivers. To track L1*L2 phase, TR first 
computes a rapidly varying (fast) model phase by scaling, according to frequency (i.e., by 
34 / 1 54 ), the C/A feedback phase obtained on the basis of a phase-locked loop tracking C/A 
phase. Because of the influence of the ionosphere, L1*L2 phase exhibits a slow drift relative to 
scaled C/A phase, with a differential rate that is usually less than 0.1 Hz for an earth-fixed 
receiver. This slow differential drift is tracked by a secondary phase-locked loop. Total model 
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phase 9m used for feedback is computed as the sum of the scaled C/A phase and a correction 
based on the output of the secondary loop. 


5.2 Correlation Function 

As in the case of known code, the correlation function for L1*L2 is obtained by applying 
an ensemble average over both noise and code to the correlation sum. Such an analysis is 
initiated by applying an ensemble average over noise (from both channels) to Eq. (5. 1). Since 
noise is uncorrelated between the LI and L2 channels, the ensemble average over the two noise 
components decouples, leading to separate averages for the two sampled voltages: 

N 

<u/(t R )> 11 = ^ < Vi(t R + ks -L x s-/s)> T ji < V* (t R + ks)>r| 2 C(9m) (5.2) 

k=-N 

where sum and averages have been interchanged. Eq. (3.13) can be substituted twice in 
Eq. (5.2), once for each averaged voltage, with appropriate labeling for channel-specific 
quantities. After the additional substitution of Eq. (4.3) in Eq. (5.2), the result can be averaged 
over code to yield 


<u,(t„)>=c,SMia£ 


a T|l G T|2 



Gi(f )G 2 (f ' ) <Ap(f ) Aj(f ')> S x df df' 


(5.3) 


where 


N 


e 2rci [cpal - q>a2 - (Pm + (f- f ')(tR + ks) + f ' T a2 - f (Xal + L x s + / s) + ^9 i] 


k=-N 


in which differenced instrumental phase is defined by 


(5.4) 


A9t=9n-9i2 (5.5) 

Again, <> denotes an average over both noise and code. Phase and delay values for LI and L2 
(i.e., with subscripts al and a2) are evaluated at the time arguments t R + k s - LxS -Is and 
tR + k s , respectively. Harmonics of the counterrotation function in Eq. (4.3) have been 
neglected under the assumption that they will average to negligible levels after the correlation 
sum. 


By comparing Eqs. (5.3) and (5.4) with Eqs. (4.7) and (4.8), one can see that an analysis 
very similar to that in the last section (5.1) can be applied to derive the following correlation 
function: 
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< u/(t R )> = A x e 2ni< P d 



Gi(f)G 2 (f) 


sin(7tfT p ) 

7C fT p 


12 


e 271 i [f Cr r - / s) +A(pi] df 


(5.6) 


where the amplitude factor is defined by 


A x — N s T p ci 


qMi qM2 sin(7t(pdT) 
1 C5V| 2 7t(p d T 


In this expression, the difference phase is defined by 

9d=(Pal(tR)-(Pa2(t R )-(pm(t R ) 


and residual delay by 


T«r ~ ^a2 (t R ) " ^al(I R ) - Lx® 


(5.7) 


(5.8) 


(5.9) 


with both time tagged at the center of the correlation interval, as indicated. To arrive at Eq. (5.6), 
use has been made of the fact that x a i and x^ are nearly equal, usually differing by less than 
0. 1 ns/s at midlatitudes at present. In addition, a phase effect of the form (p a i(/s + L x s) and a 
corresponding delay effect of the form x a i(/s + L x s), where both are caused by lagging of the LI 
samples, have been neglected under the assumption that corrections can be applied during 
processing, if necessary, to remove these accurately modeled effects. (The rates in these 
corrections can be accurately obtained from the C/A channel.) 


5.3 Discussion of the L1*L2 Correlation Function 

The L1*L2 correlation function is so similar to the code correlation function in Eq. (4.22) 
that much of the discussion in Subsection 4.4 applies here, and only the differences need to be 
outlined. Note that the filters from both channels are present in the L1*L2 correlation integral, 
rather than the filter for a single channel. This difference in filtering is important for high- 
accuracy modeling but does not change the essential character of the correlation integral. 

Because a product of channels is involved, amplitude in the L1*L2 correlation function is scaled 
downward by the RMS noise from both channels rather than only one. As discussed in 
Subsection 7.2.1, this effect greatly reduces the SNR for the L1*L2 mode of operation relative to 
the known-code mode. 

The residual delay in Eq. (5.9) is equal to the difference in the Ll-P and L2-P delays 
minus the model delay, L x , applied to the LI channel. Since this model delay must be quantized 
in steps of s, it is adjusted to the value that provides the best possible alignment under the 
circumstances, namely the alignment that places largest amplitude at the prompt (center) lag of 
the three correlator lags. Unlike the very small residual delay for code-mode processing, residual 
delay for the L1*L2 channel can be as large as half a lag (about 25 ns). 

The major component of residual phase is given theoretically by Eq. (5.8) and is equal to 
the difference of the Ll-P and L2-P phases minus the model phase. Thus, the phase being 
tracked consists of the various terms arising from Eqs. (3.15) and (3.16), but differenced between 
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channels. As in the code-mode correlation function, a small instrumental-phase component can 
emerge from the correlation integral as a result of the bandpass-phase effect, but this effect will 
be small if the receiver filters are matched between the LI and L2 channels. That is, in the P- 
codeless mode, the bandpass phase A (pi is a difference between channels. If the system filters are 
closely matched, nonlinear phase- versus-frequency variations and receiver delays tend to cancel 
between channels (i.e., A(p ( = 0 and Axi ~ 0) and therefore are significantly reduced relative to the 
case of code-mode operation. Even if the filters are exactly matched, however, the difference in 
the baseband rates for LI and L2 will prevent exact cancellation. 

For the TR receiver, one obtains the theoretical amplitude-versus-delay response shown 
in Fig. 5-l(a), assuming the dominant filter shape for both the LI and L2 channels is the 
baseband Butterworth filter described in Subsection 3.1. Fig. 5- 1(b) presents the corresponding 
phase-versus-delay response. As can be seen, phase around zero delay is nearly constant. For 
purposes of illustration, a phase offset of about 0.0053 cycles has been introduced at zero delay, 
as in the corresponding code-mode plot. Again, this offset would be removed in practice by the 
phase-locked loop and would not appear in the correlation sums. 
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Figure 5-1. L1*L2 Correlation Function for the TurboRogue Receiver 
with the Analog Front End 
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SECTION 6 


OBSERVABLE EXTRACTION 


6.1 Code Mode 

6.1.1 Phase Extraction 

A measured value of residual phase is obtained by applying an arctangent operation to the 
components of the complex correlation sum for the prompt lag ( / = 0): 


<p r = arctan[EM(uo), RE (uo)] ' 

where RE and IM respectively denote real and imaginary parts, and where "hat" over cp is used in 
Eq. (6.1) and on other observables to denote a measured value computed on the basis of 
correlation data. As implied by Eq. (4.22), measured residual phase is, on average, theoretically 
given by 


<cp r > = (pd+ <Pi 

where cpd is the difference phase in Eq. (4.10) and where (ft is the "bandpass-averaged value 
nonlinear instrumental phase not accounted for by ij. 

Measured phase at the center of the correlation interval is computed as the sum of 
measured residual phase and the model (feedback) phase at interval center: 

<P(t R ) = (Pm(t R ) + (ft 

By computing measured phase in this manner, tracking error is essentially eliminated as an error 
source in the measured delay. The resulting 20-ms phase values are fit with a quadratic function 
over successive 1-s intervals in order to reduce data volume and improve statistics. Thus, an 
aggregate phase value is extracted every second, with a time tag at the center of the 1-s interval. 
For an earth-fixed antenna, dynamics over a second are typically sufficiently restrained to allow 
a quadratic to supply an adequate fit. 

A theoretical expression for the average value of measured phase is obtained by 
combining Eqs. (3.1), (3.15), (3.16), (4.10), (6.2) and (6.3), which yields 


( 6 . 2 ) 

for 

(6.3) 


(p(t R )= <(ftt R )> = (fL-fh)tR- (phd- fL(Tg+Tt-Ts+T R ) + (p s + (p e (f L ) + 9i (6.4) 


where all quantities are time tagged at interval center, t R . As implied by this equation, the 
intentional offset in the downconversion frequency, fL-fh, which is accurately known, must be 
removed from measured phase in order to obtain a phase-delay estimate that is nearly free of 
instrumental rate offset. Measured phase delay is therefore extracted with the operation 

X(p(t R ) = [(f L -fh)t R - (p(t R )]/f L (6.5) 

Based on Eqs. (6.4) and (6.5), the theoretical estimate for measured phase delay, on average, is 
given by 
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T<p(t R ) = <X<p(t R )> = X g +X t -X S + X R + [(phd- (pi - (Ps“ (pe(fO]/fL 


( 6 . 6 ) 


where all quantities are evaluated at the center of the correlation interval and where the 
ionospheric phase, (p e , is evaluated at frequency fL ( i.e., at fu or fo). Thus, measured phase 
delay contains the following terms: geometric delay, troposphere delay, satellite-clock offset, a 
receiver-clock offset, a downconversion phase drift relative to nominal LO offset, bandpass- 
averaged receiver phase, a satellite phase offset, and ionospheric phase. 


6.1.2 Delay Extraction 

Correlation sums are accumulated over 1 s before being analyzed to extract residual 
delay. On the basis of these 1-s sums, measured residual delay can be computed from observed 
correlation amplitudes using the equation 





-k A +' A - 
' P A 0 


(6.7) 


where A., Ao, and A + are the correlation amplitudes for the early, prompt, and late lags (/ = -1, 
0, +1), respectively, and k p is a computed constant (see next paragraph) that gives residual delay 
the units of chips. The amplitude for each lag is computed on the basis of the real and imaginary 
components of the correlation sums with the projection operation: 


/N 

A; = coscpr RE[u/] + sin(p r IM[u/] 


( 6 . 8 ) 


where (p r is residual phase computed for the prompt lag according to Eq. (6.1). Based on the 
assumption that residual phase is essentially the same for all lags, this algorithm estimates 
amplitude for the early and late lags as the projection along the "direction" in the complex plane 
computed from the higher-SNR prompt lag. For the prompt lag, the algorithm yields the root- 
sum-square (RSS) of the real and imaginary components. When residual phase is large (e.g., 45 
degrees), both components of a complex correlation sum influence the computed amplitude in a 
weighted fashion, with the largest component possessing the greatest weight. Furthermore, for 
the two important side lags, the noise on those correlation sums enters the computed delay in a 
linear fashion, and a better estimate is therefore obtained at low SNR values. 

Either observed response or theoretical calculation based on an estimated system 
bandpass filter can be used to calculate the constant k p under the assumption that the residual 
delay is very small (i.e., small enough to make the amplitude difference in Eq. (6.7) 
approximately proportional to residual delay). Since the accuracy required for k p is not 
demanding, the theoretical approach has been used. Based on amplitude derivatives shown in 
Fig. 4-2(c) at the 50-ns side lags, the P-channel result is a k p value of 0.37 for the TurboRogue 
receiver. (Specifically, the theoretical value for k p is approximately 0.9/(1.28 + 1.15), where 0.9 
is the prompt amplitude obtained from Fig. 4-2(a) and the two values in the denominator are the 
derivatives at ±50 ns.) The allowed error in this constant can be relatively large as long as the 
residual delay is small. For example, an error of 10% (i.e., about 0.04) causes an error of only 
0.0001 p-chip (0.3 cm) in delay when the residual delay has a steady-state value of 0.001 chip. 
Because of the C/A-phase-driven feedback, steady-state residual delay is primarily due to the 
influence of the ionosphere and can reach 0.001 p-chip in extreme cases for an earth-fixed 
receiver, given a secondary 1-s loop of first order for tracking delay variations relative to C/A- 
derived fast feedback. 
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Residual delay values computed with Eq. (6.7) will correspond to the measure of delay 
that equalizes the early and late amplitudes. In the current version of the TR receiver with its 
unsymmetric bandpass filter in the analog front end, the resulting delay is not exactly equal to the 
delay that produces peak amplitude. For the highly digital front end described in Appendix A of 
this document, these two measures of delay would be more nearly equal due to the greater filter 
symmetry found in that implementation. (Either measure of delay could theoretically be used in 
GPS applications since they differ only by an additive constant that can be removed either by 
double differencing or by calibration.) 

Once residual delay has been extracted, 1-s measured delay is computed as the sum of 
measured residual delay and model delay at interval center: 

X(tR)=Tm(tR)+t r (6.9) 

As with measured phase, this process essentially removes any tracking error in measured delay. 
The expectation value of measured residual delay is represented theoretically by Eq. (4.24). A 
theoretical expression for the average value of measured delay is obtained by combining 
Eqs. (3.14), (4.24), and (6.9), which yields 

T(t R ) = <X> = T a s T g + t t - Xs + Tr + T e + Ti (6. 10) 

where all quantities are time tagged at the center of the 1-s interval, t R . Thus, measured delay is 
the sum of the following terms: geometric delay, troposphere delay, satellite-clock offset, 
receiver-clock offset, ionospheric delay, and receiver delay. 

To compress data volume and collect statistics, these 1-s delay values can be combined 
over longer time intervals by a least-squares fit or carrier-aided averaging. 


6.2 P-Codeless Mode 
6.2.1 Phase Extraction 

In the P-codeless mode, measured residual phase, <p r , is estimated by first adding the 
imaginary components of the correlation sums across the three lags and then normalizing to units 
of cycles: 


cp r = 



( 6 . 11 ) 


where A^ is a 3-lag-sum amplitude mapped from C/A amplitude. Unlike the code mode, 
addition across lags makes sense here since the noise on the sums is uncorrelated between lags 
(see Subappendix C.3.2). Furthermore, unlike prompt amplitude, effective amplitude for the sum 
is nearly independent of residual delay, which can be as large as 0.25 p-chip when in the P- 
codeless mode. As a result, effective SNR on the composite is about 1 .7 to 3 dB better than the 
SNR on the prompt lag alone. A small correction is needed to account for the slight dependence 
of phase on lag (about 0.015 cycles/lag, see Subsection 5.2) caused by applying lags to LI with 
its sizable baseband carrier frequency (about 308 kHz). The effect of these phase shifts can be 
accurately modeled and removed on the basis of observed LI carrier frequency, thereby 
referencing measured phase to the prompt lag. 
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As implied by Eq. (5.6), measured residual phase is, on average, theoretically given by 

<9r > = 9d+Acp I (6.12) 

where (Pd is the difference phase given by Eq. (5.8) and where A<p, is the bandpass-averaged 
value for instrumental phase, exclusive of instrumental delay effects, differenced between LI and 
L2. Measured phase for the L1*L2 channel at the center of the correlation interval is computed 
as the sum of measured residual phase and model (feedback) phase: 

9x(tR) = <Pm(tR) + 9r (6. 13) 

Again, by computing measured phase in this manner, tracking error is essentially eliminated. A 
theoretical expression for the average value of measured phase is obtained by combining 
Eqs. (5.8), (6.12), and (6.13), which yields 


<Px(Ir) = <9x> = <Pal (tR) - (pa2(tR) + ACft (6. 14) 

This equation can be reformulated through use of Eqs. (3.1), (3.15), and (3.16): 


cp x (t R ) = (Af L - Af h )t R - A(p hd - Af L (Tg+T t ) + fL2(tg 2 -T&) - fu(xL-T&) + Acps + Acpe + Acpj (6.15) 

where the superscript p refers to P channel, where all quantities are time tagged at interval center, 
and where A denotes a difference between the Ll-P and L2-P channels, in the order LI minus 
L2. As one would expect, differenced phase is the same general form as the single-channel 
phase in Eq. (6.4) but with each term differenced between channels. 

To obtain an estimate of L2 phase, the measured difference phase obtained from the 
operation in Eq. (6.13) is subtracted from the LI phase simultaneously measured in the 
C/A channel: 


CPL2 = cpu - cpx (6.16) 

where c denotes C/A channel and x denotes L1*L2 processing. Based on Eq. (6.4) as applied to 
the C/A channel and on Eq. (6.15), this measured phase, on average, is represented theoretically 
by 


(pfe = <(pL2> = (fL2 - fh2)tR - <Phd2- fl 2(Tg + T t - x| 2 + X P R2 ) + ($2+ Cp e (f L 2) + <Pl2 + V (6.17) 


where the "correction" term is defined by 

xi/ = -fu(T!,-x§i) + (cp§i-(pI 1 ) + ((pii-(pi , 2 - Acpj) (648 

All quantities are evaluated at interval center, tR, and it has been assumed that the following 
quantities are essentially identical for the Ll-C/A and the Ll-P channels: carrier frequency, 
receiver LO phase, receiver-clock error, geometric delay, tropospheric delay, and ionospheric 
phase. In practice, the satellite-clock error and satellite phase offset will be different for Ll-C/A 
and Ll-P, and these differences have been accounted for in Eq. (6.18), with a superscript to 
denote C/A or P channel. Similarly, Eq. (6.18) accounts for receiver instrumental phase. 


6-4 



Comparison of the theoretical expressions for L2 phase in the code mode, Eq. (6.4), and in the P- 
codeless mode, Eq. (6.17), shows that, under the specified assumptions, measured L2 phase 
should be identical in the two modes except for the terms in \|r. 

L2 phase delay is calculated as prescribed by Eq. (6.5), with L2 values for fL and fj,. The 
resulting phase delay, on average, is represented theoretically by 

t$(t R ) = <xJ(t R )> = Tg + X t - X P S 2 + X P R 2 + [ <Phd2- (p?2 - <P§2 - (pe(fi.2) - tfl/fn (6. 19) 

where all quantities are evaluated at interval center. Except for the instrumental terms, this 
expression is identical to Eq. (6.6). 


6.2.2 Delay Extraction 

The P-codeless technique extracts early, prompt, and late amplitudes from the correlation 
sums on the basis of the algorithm described in Subsection 6.1.2. Fig. 5-1 presents the 
theoretical dependence of amplitude on delay for L1*L2, given a 9.6-MHz Butterworth filter. 
Residual delay for L1*L2 cannot be extracted with sufficient accuracy by means of the 
algorithm, Eq. (6.7), used in the code mode because the misalignment of the Ll-P and L2-P 
correlation is not necessarily small and can be as large as 0.5 sample interval (about 25 ns). 
Instead, as illustrated in Fig. 6-1, a normalized difference of the early and late amplitudes 
(defined as the balance ratio) drives a table lookup of the residual delay. The numerical values in 
the table can be obtained through modeling and/or calibration of the amplitude response. For 
highest accuracy, separate calibration of this table is required for each receiver. The calibration 
is carried out by forcing the table lookup to produce the same measured values for P2 - PI delay 
as the code mode. To illustrate the shape of the table-lookup function, Fig. 6-2 presents 
theoretical results for residual delay as a function of the balance ratio, given the 9.6-MHz 
Butterworth filter defined in Subsection 3.1. It is estimated that the accuracy of this table lookup 
can fall between 2 and 10 cm, depending on the effort applied to calibrate the table. 

Total measured P2 - PI delay with a time tag at interval center is computed by adding the 
resulting residual delay to the feedback model delay L x : 

T'x = x r + L x s (6.20) 


A theoretical expression for the averaged value of measured difference delay is obtained by 
combining Eqs. (5.9) and (6.20), which yields 


X x — < C't x > — X a 2(tR) t a i(t R ) 


( 6 . 21 ) 


Substituting Eq. (3.14) twice, each time with appropriate subscript, yields 


t x = X§, - X§ 2 + X& 2 - X£, + X e2 - X e i + X? 2 - X?, 


( 6 . 22 ) 
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where the superscript p refers to P channel and where all quantities are evaluated at interval 
center, t R . It is assumed the geometric and tropospheric delays are identical for the two channels 
and have differenced out. In order, the remaining terms are the difference delays due to satellite 
code clocks, receiver sampler clocks, ionosphere, and instrumental delays. For accurate 
measurements of the ionospheric delay difference, calibrations are generally required to remove 
the overall delay contributed by terms other than the ionosphere. 

L2-P group delay is estimated by adding the measured L1*L2 delay to the concurrently 
measured C/A group delay: 


x p2 — ^c/a+ T'x (6.23) 


A theoretical expression for this measured delay is obtained by combining Eqs. (6.22) and (6.10) 
(with a subscript to represent LI -C/A): 


^2= < V> = ^+^-^2 + ^2 + x t 2 + ife + ?A 


(6.24) 


where the correction term is defined by 


Ta = (t!^ - T§i) - (xfi - -rfi) (6.25) 

All quantities are evaluated at interval center, t R . It has been assumed that the receiver-clock 
error and ionospheric delay are the same for LI -C/A and Ll-P. Note that the P-codeless group 
delay for L2-P in Eq. (6.24) has the same form as the code-mode measurement of the same 
quantity in Eq. (6.10), except for the differences in Eq. (6.25) due to satellite clocks and receiver 
delays. In practice, the Ll-C/A and Ll-P codes (i.e., "satellite clocks") can be misaligned by a 
few nanoseconds as they leave the satellite, and the TR instrumental delay can be on the order of 
1-ns smaller for C/A than for P. In practice, careful calibration of these two offsets should be 
carried out if one requires an accurate correspondence between the L2-P delay from the code 
mode and the L2-P delay from the P-codeless mode. 
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SECTION 7 


SYSTEM-NOISE ANALYSIS 


7.1 Code Mode 

7.1.1 Signal-to-Noise Ratio 

The signal-to-noise-ratio analysis presented in this section is based on the analog front 
end described in Section 3, with two-level sampling. The voltage SNR for a complex correlation 
sum is defined as the maximum correlation amplitude (i.e., RSS of components for prompt lag) 
divided by the RMS noise on either component. For perfect lock in phase and delay in which 
residual phase, phase rate, and delay are all equal to zero, Eq. (4.22) predicts that the maximum 
amplitude Am at the prompt lag is given by 


Am=N s T pCi ^ 


r G(f) 

sin(7tfT p ) 

L <0 

TtfTp 


e 27ti&p, df 


(7.1) 


where 8(pj is the nonlinear phase-frequency response of the receiver filter but with an integral - 
averaged value subtracted (by the DPLL) so that the integral is real. 

The integral, which is a measure of the maximum correlation of the filtered baseband 
signal with the pure P-code sequence generated by the receiver, is approximately proportional to 
the square root of signal power at baseband. Based on this proportionality, one can set up the 
following relationship between the correlation integral and signal power: 


Ip 

Grj 
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I CXf } 

sin(7tfT p )' 

1 ) 

* -oo 

tt fT p 


2 


e 2iti&p, df 



(7.2) 


where is proportional to the average noise power Pn entering the sampler after receiver 
filtering, and where P s is the signal power per quadrature component that would have entered the 
sampler if the receiver had not filtered out any frequency components of the signal. The factor of 
2 under the square root accounts for the two quadrature components, while the factor Of 
accounts for loss in signal power due to receiver filtering. Numerical calculations indicate that 
ctf is approximately 0.9 for the P channels and 0.99 for the C/A channel, given the analog front 
end described in Section 3. This method of representing the maximum correlation has been 
adopted in order to explicitly formulate Eq. (7.1) in terms of a measure of signal power that is 
independent of the shape of TR filters and to explicitly account for the loss due to receiver 
filtering. 

When Eq. (7.2) is substituted in Eq. (7.1), the maximum amplitude becomes 


V 


p CtfPs 

Pn 


Am = Ns Ci qM 


(7.3) 



The ratio, 0Cf P s /Pn, can be thought of as the single-sample power SNR per quadrature component 
at baseband after accounting for filtering but neglecting quantization loss. 

Maximum voltage SNR is computed by dividing maximum amplitude by the RMS 
noise, G u , on the either component of the complex sum: 

SNR V = ^ M - (7.4) 

<*u 

Based on Eqs. (7.3) and (C.25), the maximum voltage SNR in the case of quadrature sampling 
becomes 


SNR v = L c L s y2N^S 


(7.5) 


where N s is the number of samples correlated, not counting samples lost during the dead time 
between correlation intervals, and where the factors L c and L s are the losses due to quantization 
of counterrotation sinusoids and sampling quantization, respectively. The loss caused by three- 
level quantization of the counterrotation sinusoids is defined as 

L c = -pL- = = 0.96 (7.6) 

iL5 fL5 


or 0.35 dB. The sampling quantization loss is defined by 



(7.7) 


Based on Eqs. (B.15) and (B.16), this loss is 

L s = yT = 0.798 (7-8) 

or 2 dB for the two-level sampling implemented in the current TR receiver. 

The term N s in Eq. (7.5) can be thought of as the SNR increase due to N s independent 
samples. Ideal sampling is assumed; i.e., the sample rate and bandpass are such that the noise on 
sample points is essentially independent between points. For nonquadrature sampling, the result 
in Eq. (7.5 ) decrea ses by V2. Total proce ssi ng loss due to the receiver is approximately equal to 
L c L s Vctf V18/20 or 0.96 x 0.798 xV0.9 x V 18/20 = 0.69 (3.2 dB), where 18/20 accounts for 
samples lost during the dead time. 

A nominal SNR for the TurboRogue receiver in the Ll-P channel is calculated as 
follows: The total noise power at RF at insertion (i.e., including antenna gain but before front- 
end amplification) is given by 

PN=k B T s W (7.9) 

where ks is Boltzmann's constant, T s is system (total) noise temperature, and W is the effective 
double-sided noise bandwidth. For a nominal system temperature of 200 K and a noise 
bandwidth of 19.2 MHz, the noise power becomes 
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Pn= 1.38 x 10' 23 w/K-Hz x 200 K x 19.2 MHz 


(7.10) 


ex 

Pn= 5.3 xlO' 14 watts (7.11) 

The noise bandwidth of 2*9.6 MHz = 19.2 MHz reflects the ultimate double-sided passband of 
the front end. 

According to GPS specifications, nominal signal power collected by an antenna with 
0-dB gain is -163 dBW for the Ll-P signal so that the signal power at insertion at RF for a 
3-dB-gain antenna is 


Ps = 10' 16 watts (7.12) 

Since the receiver amplifies signal and noise by the same factor in transit from RF to baseband, 
and since downconversion does not change the ratio of signal to noise, the baseband power ratio 
in Eq. (7.5) is the ratio of the numbers from Eqs. (7.1 1) and (7.12): 

£^ = 0.0019 (7.13) 

For a 1-s integration time consisting of 50 intervals of 18-s duration and a sample rate of about 
20 MS/s, Eq. (7.5) predicts a nominal TurboRogue voltage SNR of about 190, given the losses 
specified above. For the latest version of TurboRogue, typical zenith SNRs for Ll-P are about 
600, which is larger due to the effect of a 7-dB antenna gain at zenith, 1 10-K system temperature, 
and higher-than-nominal signal strength. (Signal strengths are reportedly 3 to 5 dB greater than 
nominal.) 

Similar calculations can be carried out for the C/A and P2 channels. As a consequence of 
differences in transmitted power, the nominal C/A SNR is approximately V2 higher and the 
nominal L2-P SNR approximately V2 lower than the Ll-P result. 


7.1.2 System-Noise Error in Phase 

The noise on the correlation sums propagates to noise in extracted phase. By propagating 
noise through the arctangent operation in Eq. (6.1), one can show that the system-noise error in 
phase (la, in cycles) is given by 


~ 27t SNR 20 


(7.14) 


where SNR 2 o , as computed in Eq. (7.5), is the 18-ms voltage SNR entering the arctangent. For a 
1-s observable, the phase error can be computed by substituting the SNR from Eq. (7.5) in 
Eq. (7.14) but increasing the error by an additional factor of 1.5. The factor of 1.5 is the 
degradation suffered relative to a straight average when a quadratic fit is applied to fifty 20-ms 
phase values extracted by the phase-locked loop over the second. Based on the nominal 1-s SNR 
of 190 presented in the Subsection 7.1.1, the system-noise error in Ll-P phase therefore becomes 
0.00126 mcycles (0.024 cm) for a 1-s integration time. 
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After further increasing the integration time by means of an additional quadratic fit over 
Nt seconds (Nt > 2), the system-noise error is decreased by a factor of / 1.5. The 
additional factor of 1.5 is the degradation due to the second quadratic fit. 


7.1.3 System - Noise Error in Group Delay 

System-noise error in delay is calculated by propagating errors through Eqs. (6.7) and 
(6.8). To first order, one can easily show onthe basis of Eq. (6.8) that the amplitude errors for 
the side lags are independent of the error in <p r , with RMS values given by 

Ga = (sincp r ) 2 Gu + (cos(p r ) 2 G^ (7.15) 


which leads to 


cJa = Ou (7.16) 

where G u is the 1-G error in either component of the complex correlation sum. This result 
assumes negligible correlation between the noise on the quadrature components of the correlation 
sums as derived in Eq. (C.23). 

Noise on the three amplitudes propagates through Eq. (6.7) to create a delay error given 
by 


(7.1 

Ao Ao 

where 8 denotes deviation from the actual value. Because residual delay, t r , is small during in- 
lock code-mode operation, the 5Ao term in this expression can be neglected. The RMS of the 
delay error in Eq. (7. 17) becomes 


G T = V2k p ^[l-A(2)] 1/2 


(7.18) 


where A is the correlation of noise between lags with / - 1"= 2 (see Eq. (C-26)). If the prompt lag 
has peak amplitude, Eqs. (7.4), (7.16), and (7.18) can be used to show 



(7.19) 


again with units of chips. The correlation, A(2) has been set to zero, as explained in 
Subsection C.2.2. For the nominal 1-s SNR of 190 obtained above for Ll-P and kp = 0.37, the 
delay error is equal to approximately 0.0028 p-chip (8.3 cm) for the Ll-P channel for a 1-s 
integration time. For the aforementioned observed zenith SNR of 600, the delay error is 2.6 cm. 

If delay is averaged over a longer time inter val o f length Nt seconds by means of carrier- 
aided averaging, the error decreases by a factor of VNt- For 5 min, the error is a factor of V300 
smaller or 0. 15 cm for the maximum observed SNR. 
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7.2 P-Codeless Mode 
7.2.1 Signal-to-Noise Ratio 

The definition of signal-to-noise ratio for the L1*L2 channel is the same as that for the 
code channels presented in Subsection 7.1.1. By using an analysis that parallels the steps in that 
subsection, one can derive a corresponding expression for L1*L2 voltage SNR. Based on 
Eqs. (5.6) and (5.7), maximum correlation amplitude is given by 


Am — N s Tp Cj 


qMi qM2 

°r| 1 ^T|2 


Gi(f)G 2 (f) 


’sin(TifTp)" 

JtfT„ 


e 2jti8(pi <Jf 


(7.20) 


Analogous to Eq. (7.2), the correlation integral is related to satellite power by 
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(7.21) 


r% 

where is proportional to the average noise power Pnk entering the sampler for channel Lk 
after receiver filtering, where P s k is the average signal power at baseband for the P signal in 
channel Lk, as defined following Eq. (7.3). The phase, 8tpi, corresponds to the definition for 
Eq. (7.1) but is based here on the difference of receiver phase between P channels. Again, oif 
represents the loss due to receiver filtering. Numerical evaluation based on the adopted 
Butterworth filter described in Subsection 3.1 indicates that 0Cf is approximately equal to 0.9. 

When Eq. (7.21) is substituted in Eq. (7.20), the maximum amplitude becomes 

Am = N s ci q M i qM2 V ^™ (7.22) 

Note that Eq. (7.22) has a factor of 4 under the square root rather than the factor of 2 in Eq. (7.3). 
This increase accounts for the fact that four quadrature components contribute to the L1*L2 
correlation rather than two. 

Maximum voltage SNR found at perfect alignment of the LI and L2 signals is obtained 
by substituting Eqs. (7.22) and (C.35) in Eq. (7.4), which yields 


SNR X = L c 


(7.23) 


where all quantities have the same definitions as in Subsection 7.1.1, and the subscript numbers 
refer to the Ll-P and L2-P channels. Since the counterrotation sinusoids are again quantized to 
three levels, the quantization loss, L c , again has the value derived in Eq. (7.6). The sampling 
quantization loss, however, is given by 



qMi qM2 


(7.24) 
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Based on the results in Eqs. (B.15) and (B.16), this loss becomes 

L s = — = 0.637 

Jt 

given two-level sampling quantization in both the LI and L2 channels. 


(7.25) 


The nominal L1*L2 SNR for the TurboRogue receiver with the analog front end can be 
calculated on the basis of the parameter values in Subsection 7.1.1. When the power ratio for PI 
in Eq. (7.13) and the above losses are substituted in Eq. (7.23), the maximum voltage SNR on the 
prompt lag for a 1-s integration is approximately 4.7 where the value for Ps 2 /Pn 2 has been set 
equal to 1/2 of the same quantity for LI to reflect the nominal decrease in transmitted power for 
L2. Three-lag combining improves the effective SNR by 1.8 dB for phase (see Subsection 7.2.2), 
which changes the above SNR value to 5.8 for phase tracking and estimation. 

Actual SNRs are considerably better than these nominal SNRs. With the current version 
of TR, the peak observed 1-s SNRs are approximately 75 instead of 4.7, and 90 instead of 5.8. 

As mentioned in Subsection 7.1.1, the increase is a result of lower system temperature, higher 
antenna gain, and larger-than-nominal source strength than assumed in the nominal calculation. 

A useful ratio of SNRs is defined by 


R _ SNRvi SNRv 2 
cx SNR X 


(7.26) 


where SNR vie is the 1-s voltage SNR at the prompt lag for channel Lk in the code mode and 
SNR X is the 1-s voltage SNR at the prompt lag in the P-codeless mode. This ratio, which is a 
good check of receiver health, should be nearly constant over a satellite pass. Eqs. (7.5) and 
(7.23) lead to the expression 


Rex — Lc V 2 otfN s 


which becomes 


R cx = 0.96 V 2 x 0.9 x 20.456 x 10 6 MS/s x 0.018 s x 50 - 5500 (7 ' 28) 

for the TurboRogue receiver. It has been assumed that a 1-s integration contains fifty 18-ms 
intervals. Observed ratios will deviate slightly from this approximate value. Unlike the 
assumption behind this result, noise is slightly correlated between 20.456-MHz sample points. 
Further, the observed ratio is based on an estimate of peak prompt L1*L2 amplitude, which is 
obtained by combining correlation amplitudes for the three lags and is therefore corrupted by 
errors in modeling amplitude dependence on lag. (In practice, the amplitude at the prompt lag in 
the LI *L2 channel is not peak amplitude; peak amplitude must be inferred from the three 
amplitudes.) Preliminary measurements have produced an observed ratio of approximately 5000 
(C.E. Dunn, 1994). 
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7.2.2 System-Noise Error in Phase 

As discussed in Subsection 6.2.1, residual phase is extracted by means of a sine phase 
extractor applied to the sum of the Q correlation sums for the three lags. Relative to peak 
amplitude, the combined amplitude is approximately 1.0 + 2*0.563 = 2.13, given side 
amplitudes of about 0.563, as suggested in Fig. 5-1. Since noise on the correlation sums is 
uncorrelated between lags (see Subappendix C.3.2 of this document), noise on the combined 
sums is i3 larger than the noise for any given lag. Thus, combining lags improves voltage SNR 
by about 2.13/1.732 = 1.23, or 1.8 dB, relative to the phase error at the peak amplitude alone. 
(Since the prompt lag is not necessarily at peak amplitude, combining lags generally improves 
SNR by even more than 1.8 dB relative to the prompt-lag SNR — by as much as another 1.2 dB 
for a total gain as large as 3 dB.) To compute system-noise error in phase, voltage SNR from 
Eq. (7.23) can be substituted in Eq. (7.14) and adjusted by the 1.8-dB SNR increase. Based on 
the nominal 1-s SNR of 4.7 computed at the correlation peak amplitude in Subsection 7.2.1, the 
1-s system-noise error in L1*L2 phase becomes 0.028 cycle or 0.67 cm for a 1-s integration time, 
after the aforementioned 1.8-dB increase in SNR. (There is no loss of 1.5 for L1*L2 phase since 
C/A-aided averaging compresses L1*L2 phase rather than a quadratic fit.) Again, in practice, 
actual errors are considerably smaller than this nominal error. For the zenith SNR of 90 quoted 
in Subsection 7.2.1 for phase, the 1-s error is approximately 0.04 cm. 


7.2.3 System-Noise Error in Delay 

As discussed in Subsection 6.2.2 and illustrated in Fig. 6-1, delay extraction is based on a 
prompt-amplitude-normalized difference of early and late amplitudes. For the purpose of 
estimating system-noise error in delay, the table lookup in Fig. 6-1 can be approximated by 
Eq. (6.7) with a k p value of 0.37. When the system-noise errors in amplitude are propagated 
through Eq. (6.7), one obtains the error presented in Eq. (7.17). The RMS (i.e., l-o) error is 
obtained by squaring the error in Eq. (7.17), averaging over noise, and taking the square root. 
Since the amplitude errors due to the system noise are uncorrelated between lags, the 1-a delay 
error becomes 



a ^ kp Ao 
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(7.29) 


in units of p-chips. 


This expression for delay error is a weak function of the imbalance between early and late 
amplitudes. An upper limit for a given SNR is obtained by letting A+ = Ao- Based on the 
amplitude-versus-delay plot in Fig. 5-1, the delay error in this case becomes 


G x <k p 


1.15 

SNR X 


V2 + 0.51 


(7.30) 


in units of p-chips, where 1.15 accounts for the decrease of the prompt amplitude relative to peak 
amplitude and 0.51 is approximately the maximum value for the normalized difference of the 
early and late amplitudes. Eqs. (7.4) and (7.16) have been used to convert to SNR X , the peak 1-s 
SNR. Thus, an upper limit on delay error is given by 




0.67 

SNR X 


(7.31) 
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in units of p-chips. The approximate k p value of 0.37 has been assumed. 

When the averaging interval has a length of Nt seconds, the error is reduced by a factor 
of VNt . For a 5-min integration and the nominal 1-s SNR of 4.7 computed in Subsection 7.2.1, 
this upper limit for system-noise error in delay becomes 

a < 0.67 1 

T 4.7 V300 (7.32) 


or 


a T < 0.0082 p-chip (25 cm) (7.33) 

The corresponding lower limit on delay error (i.e., when early amplitude equals late amplitude) is 
about 29% smaller than this value. When computed on the basis of the aforementioned observed 
maximum delay SNR of 75, this nominal delay error reduces to 1.6 cm. 
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APPENDIX A 


EXAMPLE OF A HIGHLY DIGITAL FRONT END 

The front end presented in the text can be improved through greater use of digital 
processing to offer greater stability, lower cost, improved delay accuracy, smaller size, and less 
weight. As an example of a possible future improvement relative to the analog front end, 
therefore, a highly digital front end has been designed for the TurboRogue receiver. This 
appendix functionally and analytically describes for that particular front-end design the 
processing that takes the signal from reception at RF down to its sampled form at baseband. The 
proposed front end is based on downconversion sampling, a technique in which the signal is 
sampled at RF with a rate commensurate with a selected midband RF frequency, thereby causing 
that particular frequency component to be aliased to zero frequency. The spectrum surrounding 
the selected component is aliased to baseband as though it had been downconverted with a 
mixing signal having the selected frequency. In this manner, sampling and downconversion 
from RF are carried out in one step, avoiding the mixer and filter associated with standard 
downconversion. In the design proposed here, the effective LO frequencies and the final 
baseband sample rates are the same as those in the analog front end described in the text. 


A.1 Functional Description of a Highly Digital Front End 

A high-level functional block diagram of the example front end is shown in Fig. A-l. 

The RF signal collected by the antenna is passed through an "anti-RFI" filter that is sufficiently 
wide to pass both the LI and the L2 bands. After broadband amplification, the signal is power- 
divided into LI and L2 branches, and each branch is passed through a filter (defined as channel 
filter) centered at the appropriate L-band carrier. Width of these filters depends on the 
application, with narrow widths selected to provide good RFI rejection and wider widths to allow 
better delay calibration. (The narrower widths will also provide slightly better SNR performance 
because of the "oversampling effect" (Beaulieu, 1988) in the case of 1-bit sampling.) Figure A-l 
is an example of a wide filter where the width is set equal to the largest advisable value, namely, 
the sample rate for each quadrature channel. After filtering, each channel is sign sampled (1-bit 
quantization) at RF at a very high rate (900.064 MHz for LI and about 981.888 MHz for L2) 

The sample points are passed to an alternate-sample demultiplexer (demux) that separates 
the C (defined as cosine) samples from the S (defined as sine) samples, thereby producing two 
data streams, each at half the original rate. A sign correction is applied to alternate samples on 
each stream to remove the sign flips described below. A sum-and-dump operation (22 points for 
LI and 24 points for L2) is then applied to reduce the data rate to 20.456 MHz per quadrature 
component for both the LI channel and the L2 channel. As a result of the sum-and-dump 
operation, the number of bits required to represent the samples is increased from 1 to 6. For 
example, sum-output values can range between -22 and +22 for the LI channel, in principle. 
Since the TR baseband processor is designed to accept only two- or three-level samples, the 6-bit 
samples are requantized to three levels: (-1,0, +1). This requantization requires specification of 
two integer discriminator values that set the decision levels for transforming each sum value to 
-1, 0, or +1 . These discriminator values are set to values that minimize SNR loss and are 
consequently a function of channel bandwidth. When the largest advisable channel bandwidth is 
selected (e.g., 450 MHz for L2), the SNR loss due. to requantization is minimized when the 
discriminator levels are set so that X < - 3 => - 1 ; |x| < 3 => 0 and X > 3 => + 1. 
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Figure A-l. Block Diagram for the Example of a Highly Digital Front End 












A.2 Explanation of the Highly Digital Front End 


As mentioned above, the two RF sample rates are set to values nearly commensurate with 
their respective carriers so that each carrier will alias to nearly zero frequency as a result of 
sampling. Thus, the sampling process also is a downconversion process. In addition to being 
nearly commensurate with the carrier, the sample rate is selected so that the sample points 
alternate between quadrature components of the signal. 

These features of sampling downconversion are schematically illustrated in Fig. A-2 for 
the L2 channel. In this example, the sample clock has a rate of F s = 4/5 fh, where fh is the desired 
effective downconversion frequency. (In TR, the effective downconversion frequency is close to 
the carrier frequency, with an offset of 308 kHz for LI and 240 kHz for L2, as in the analog front 
end.) Figure A-2 shows the location of several sample points on a sinusoid with a frequency of 
fh. For that particular frequency, the alternate points labeled "C" for cosine all have the same 
amplitude value except for an alternating sign flip. Thus, after a correction for the sign flips, the 
C samples for the fh component are equal to a constant or, as commonly described, alias to zero 
frequency. The same observations can be made for the indicated alternate sine ("S" ) points. 
Figure A-2 also illustrates the quadrature relationship of alternate points. Note that adjacent C 
and S points differ by a quarter of a cycle in the fractional cycle. This relative phase of a quarter 
cycle reveals the quadrature relationship between pairs of adjacent sample points, as 
mathematically shown in the next subsection (A.3). In this manner, RF sampling produces the 
two quadrature channels as interleaved samples. 

A similar graphical analysis can be applied to a frequency component displaced from fh- 
For the (zero-Doppler) carrier frequency fL, for example, such an analysis would show the 
sample points shifting in phase relative the cycle boundaries in a way that makes the frequency 
based on observed sample values appear to be fL-fh- Thus, all frequency components in the 
sampled band will effectively be lowered in value by fh- A schematic illustration of the 
frequency-domain consequences of sampling downconversion is shown in Fig. A-3. (This figure 
neglects the distortion inflicted on the noise spectrum by the subsequent two-level amplitude 
quantization described below.) Note that this method results in a so-called "double-sideband" 
downconversion, in which half of the baseband spectrum has negative frequency. These 
characteristics of sampling downconversion are demonstrated mathematically in the next 
subsection. 

As implied by Fig. A-l, the frequency subsystem has to generate only three frequencies, 
the ultimate baseband sample rate (commonly referred to as sample clock) of f s = 20.456 MHz 
and two RF sample clocks obtained by simple multiplication of f s : F s = 44f s for LI and 48f s for 
L2. As in the analog front end, f s is offset from 20.46 MHz by 4 kHz in order to ensure effective 
incommensurability with the P chip rate. Also paralleling the analog implementation, the 4-kHz 
offset produces an offset in both of the effective downconversion frequencies (relative to the 
zero-Doppler carrier). Again, the three required offsets can all be changed by one adjustment in 
f s , without changing any multipliers. 

The example implementation shown in Fig. A-l assigns each channel filter the largest 
advisable width, namely a width approximately equal to the RF sample rate per quadrature 
channel. In selecting these filter widths, one finds that there is a trade-off between RFI rejection 
and receiver-induced errors in phase and delay. If the LI filter width is reduced from 450 MHz 
to 50 MHz, for example, then the effects of RFI outside that 50-MHz band will be greatly 
reduced. (In addition, subsequent oversampling at a 450-MHz rate would help to recover some 
of the SNR loss caused by two-level quantization.) Accuracy, on the other hand, would be 
degraded, since analog filters introduce errors and instabilities in phase and delay that increase as 
bandwidth decreases. Such errors can be reduced by stable designs and calibration, but they can 
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Figure A-3. Schematic Illustration of Sampling Downconversion in the Frequency Domain 





be more reliably reduced by increasing the bandwidth. For example, the delay through a 
450-MHz Butterworth filter is a few nanoseconds, depending on the number of poles. Given the 
reasonable goal of 5 % calibration and stability, one could reduce the channel-filter delay error to 
the order of 0.2 ns. In contrast, the delay error would be about nine times larger for a 50-MHz 
filter. 


A.3 Mathematical Model for Sampling Downconversion 

This subsection presents a mathematical formulation of the sampling downconversion 
process outlined in the previous subsections. After application of the antenna filter and the 
complex channel filter, G c , to the received signal in Eq. (2.19), the RF voltage becomes 


V F (t) = e 2lci (fLt-f L t) Gc(v d )G A (v d )Gs(f)A p (f )e^ i [f(t-T)+<Ps+«Pel df + c.c. + r| F (A.l) 


where T|f is the RF noise passed by the channel filter and where the argument, Vd, given by 
Eq. (3.3), is the Doppler-shifted frequency component at RF corresponding to f. The carrier 
factor has been pulled outside of the integral so that it can be reformulated in a manner that 
anticipates downconversion: 


V F (t)= e 2lci lfht+(fL-fh)t-fLT] Gc(Vd)G A Gs(f)Ap(f)e 2lti [f(t-'c)+<ps+(pddf + c.c.+ % (A.2) 


where fh is the anticipated effective downconversion frequency for either LI or L2 as discussed 
in the preceding subsection and as defined in Fig. A-l. 

The next step in processing is to sign-sample each of the two channels at RF. Following 
the analysis in Subsection 3.3 for sign sampling in the analog front end, one obtains the 
following expression for a sample after an ensemble average over noise: 


<V F (t Rk )> T1 ~ S!le 27ti [fhtRk+(fL-fh)tRk-fLi:'] 


[ 


G c (Vd) G A G s (f ) A p (f )e 27ti [fte- r')+ cps + (Pel df 


+ c.c. (A.3) 

where tRk represents the kth sample time according to the RF sampler clock. The approximately 
equal sign denotes the neglect of the higher order terms produced by sign sampling, as discussed 
in Appendix B of this document. For the example design in Fig. A-l, omission of the higher 
order terms is an excellent approximation since the wide noise bandwidth results in a very small 
single-sample SNR (e.g.. S/a < 0.04). Again, delay has been redefined, as in Eq. (3.7), to 
incorporate the sampler clock error that enters when true time is replaced with sampler time. 

Only the L2 samples will be analyzed in detail since the LI analysis follows an identical 
course. In the proposed implementation, the RF sample rate for L2 is 4/5 fh. The sample rate 
becomes 981.888 MHz since fh = 60*20.456 MHz. The sample times can be represented in 
terms of this sample rate by 



tRk = k , k= 1,2,3, - 

41h 


(A.4) 


These sample times will be separated into alternate points: 


t c R„ = 2n5J- 

4f h 


and 


t s R „ = (2n-l)5.f 
4 th 


(A.5) 


(A.6) 


in which n = 1,2, 3..., and where the superscript C denotes cosine (real component) and the 
superscript S denotes sine (imaginary component). 

For the C samples, tRk is replaced with tL in Eq. (A.3) to yield 


<V F (tL)> n «^. e 2itit2.5n+(fL-fh)t| n -fLf]j[ Gc(Vd) G A Gs (f ) A p (f ) e 2jt ' [ f (tRn- ) + <Ps+ <p e ] df 


+ C.C. 

which can be rewritten as 


(A.7) 


<Vp(tS n )> T1 *(-l) n ^ 

Or| 


e 2jti [(f L - fh)tL-fLt] 


l 


Gc(Vd) G A Gs(f )Ap(f ) e 2rci[f(tRn-T’) + <Ps+ 9el df 


+ c.c. (A.8) 

where %' and (j) e are evaluated at time tR n . Thus, except for a sign flip on alternate samples, the 
carrier of the C samples is effectively downcon verted by frequency fh, as indicated by the 
disappearance of the fj,t term in phase. After the receiver removes the sign flip, the C samples 
can be rewritten as 


<V F (tL)> 11 



Re[G c (Vd) G A G s (f )Ap(f)e2^i[(fL-fh)t c R„-fLT'+f(t c Rn -f )+cps + (p e ]] df 


(A.9) 


where Re[ ] denotes real component. 

In a similar fashion, one can show that the S samples are given by 
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<V F (tL)>n «(-l) n ^e 2jti f( f L-fhXL-fL^'-1.25] I 

^ L 


Gc(V d )G A G s (f )A p (f ) e 27ti[f(tL-T , )+(ps+ <p e ] df 


+ c.c. (A. 10) 

where t' and (|)e are evaluated at time t^ n . Since sign flips, (-l) n , are removed in processing and 
since 

e 2jii(-1.25) = _i (A. 11) 


Eq. (A. 10) becomes 


<V F (tL)> n 



Im[Gc(v d )GAGs(f)A p (f)e 2ni [(fL-fh)tL-fLT' + f(tL-t')+(ps + cpe]]df (A. 12) 


where Im[ ] denotes imaginary part. 

Except for the small separation of the C and S sample times, Eqs. (A.9) and (A. 12) are 
the real and imaginary parts of the same quantity. This time separation, At F , is equal in 
magnitude to the sample spacing at the initial RF sample rate (f s = 4 fh / 5 = 48 * 20.456 MHz 
for L2) and is close to a nanosecond (about 1.0184 ns for L2) for the proposed implementation. 
As shown in the next two paragraphs, the effect of the time difference is negligible in carrier 
phase and is known and removable in group delay. Thus, when the samples represented by Eqs. 
(A.9) and (A. 12) are combined as a complex number, they can be regarded as a complex sample 
of the baseband signal, with a sample rate of 24*20.456 MHz per quadrature component for L2. 

The effects of the time separation of the C and S components can be modeled as follows: 
The carrier phase terms in the S-sample expression in Eq. (A. 12) are defined by 


tPc(tRn) = (fL-fh)tRn “ flAf tL) + (Ps(tRn) + (Pe(tRn) 


(A. 13) 


and can be rewritten in terms of the C-sample times as 


Cp c (tL) « Cp c (tR„) - A(p F 


(A. 14) 


where the S "phase error" is given by 


Acp F s (ft - fh - fiA' + (p s + (p e ) At F (A. 1 5) 

Since fh is offset by 240 kHz from the carrier for the L2 channel, since the maximum Doppler 
shift for an Earth-orbiting satellite is on the order of 50 kHz, and since the maximum ionosphere 
rate is on the order of 1 Hz at L band, this phase error for L2 is bounded by 


Acp F < (240 kHz ± 50 kHz± 1 Hz) x Ins 


(A. 16) 
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which gives the upper limit 


|a<|)f| < 3 x 10‘ 4 cycles (0.07 mm) (A. 17) 

Note that this phase error is nearly constant, with a time-varying component on the order of 
50 kHz*l ns = 50 pcyc = 0.012 mm or less. A similar negligible error characterizes the LI 
channel. Thus, the C and S samples are essentially simultaneous with respect to phase. 

The terms in the phase of Eq. (A. 12) that can potentially effect measured group delay can 
be rewritten in terms of C-sample times as 


tftn - T'( tRn) ~ tRn “ t'( tRn) - ATf 


(A. 18) 


where the S "delay error" is given by 


At f = (1 - i') At (A. 19) 

Since maximum delay rate is on the order of 40 (is/s for the extreme case of an Earth-orbiting 
receiver, the contribution of the delay-rate term to this error is negligible (on the order of 
4 x 10~ 14 seconds). Thus, the S delay error is a constant bias equal to the sample spacing (~1 ns). 
This is the delay error that would be present in the S data if one processed the S samples as 
though they had been recorded at the C-sample times. As explained below, the C and S sample 
pair is processed as though both components were recorded at the C sample time. The single 
delay value that is extracted from this combined data will be an average of the C-sample delay 
and the biased S-sample delay. Consequently, the measured delay will be offset from "true 
delay" by a constant bias equal to Atf/2 (0.5555 ns for LI and 0.5092 ns for L2). In 
applications involving extremely accurate clock synchronization, this delay bias might be 
important, but it can be accurately corrected, if necessary. In applications using double 
differencing to remove receiver "clock" effects, this error differences out and is of no 
consequence. 

The above analysis shows that C and S sample pairs can be processed as simultaneous 
quadrature samples at time t^„ , provided a small, accurately known delay correction is applied 
to measured group delay after the data is processed. Thus, in the analysis that follows, adjacent 
C and S samples will be modeled as simultaneous quadrature samples. A composite expression 
for the complex sampled signal is obtained by replacing t| n by t^ in Eq. (A. 12) and forming a 
complex quantity with the real part occupied by Eq. (A.9) and the imaginary part occupied by the 
adjusted version of Eq. (A. 12). Such a combination yields 


<Vp(tW>,=^ 


G c (Vd) G a G s (f ) A p (f )e 27ti [(fL+f)(tRn-r')+(ps + <Pe - <Phl df 


(A.20) 


where the factor of two has been absorbed in the receiver filter term Gc- The downcon version 
phase containing fj, has been placed in <Ph , as defined in Eq. (3.1), but evaluated at time tp> n with 
(phd set equal to zero. 

By comparing Eq. (A.20) and Eq. (3.6), one sees that sampling downconversion can 
produce a sampled signal that has the same basic form as a signal subjected to the analog 
downconversion, filtering, and sampling described in Subsection 3.3. In the proposed sampling 



downconversion approach, however, the noise bandwidth and the sample rate are much larger 
(e.g., 490.944 MHz per component for L2 versus 20.456 MHz). The next subsection models the 
sum-and-dump operation used to reduce the high initial data rate to 20.456 MHz. 


A.4 Mathematical Model for the Sum-and-Dump Operation 

The sum-and-dump step compresses successive groups of Nr points (22 for LI and 24 for 
L2) to one point using the operation 


Nr-1 

V(t R ) = £ V F (t&0 (A.21) 

k = 0 

where t R is the center of the sum interval. The average value of this sum is found by taking an 
ensemble average over noise: 


n r -i 

<V(t R )> n = X <V F (t& k )> l , 

k = 0 

Next, represent time relative to interval center, t R , by 

tk=t R + (k- Nr/ 2+ l/2)s' 

where s’ is the RF sample separation given by s' = 2/F s (approximately 2 ns). Now expand the 
time-dependent phase quantities in Eq. (A.20) about center time to obtain 

(f L +fX tRk-T') + (pS + <pb-(ph = 

[fL+f][t R -T’(t R )] + cps(t R ) + cpb(t R ) -<ph(t R ) + f d |k-N R /2+ 1/2] s' (A.24) 

where f d is the Doppler-shifted baseband frequency defined in Eq. (3.4). When Eq. (A.20) and 
(A.24) are substituted in Eq. (A.22) and the sum and integral interchanged, one can show 


(A.22) 


(A.23) 


<V(W>, = Nr ^[ 

J-oo 


G R (fd) G A Gs(f)Ap(f ) e 2rci[( f L+f)(t R -T )+ IPs + <Pe-<Ph] df 


(A.25) 


where all phase terms are now evaluated at time t R and where the receiver filter is defined for this 
implementation by 


G R (f d ) = G c (v d ) 


sin[N R 7tfdS'] 

N R sin[7tfds'] 


(A.26) 


The relation 
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sin[N R 7tfs] 

sin[7cfs] 


(A.27) 


has been used to evaluate the sum over time. 

The sum over Nr samples increases the range of sample values from two possible values 
(-1 and +1) to - Nr to + Nr, thereby increasing the number of bits required to represent a sample 
from 1 to 6. Since the TR digital signal processor will only accept two- or three-level 
quantization, the sum output is requantized to three levels (-1, 0, +1), as described in 
Subappendix B.1.4. In the low-SNR limit that applies to TR processing, this requantization does 
not change the ensemble-averaged form for the signal found in Eq. (A.25), except for a scale 
factor. Three-level requantization causes a loss in voltage SNR of about 10% (0.9 dB). 

Note that the theoretical expression for the ensemble-averaged baseband samples is the 
same form for the analog front end modeled by Eq. (3.6) and the digital front end modeled by 
Eq. (A.25) except for a scale factor. For this reason, the analysis in the text applies directly to 
digital-front-end samples. 
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APPENDIX B 


QUANTIZATION EFFECTS 


B.l Sampling Quantization 

Amplitude quantization carried out during the analog-to-digital conversion of the signal 
can introduce distortions and decrease SNR. Two-level and three-level quantization are analyzed 
in this appendix. 


B.1.1 Basic Principles 

The complex analog voltage at baseband in Eq. (3.2) can be expressed in terms of its real 
and imaginary parts by 


V = V x +iV y 


(B.l) 


and in terms of signal and noise by 


V = (S x +T| x ) + i (Sy+ri y ) (B.2) 

where S and r| represent the signal component and noise components, respectively, with x and y 
denoting real and imaginary parts, respectively. The following analysis applies to both real or 
imaginary parts but will be labeled x. 

In its ideal form, sampler quantization of voltage can be modeled by a quantization 
function Q(V X ) that is a piecewise-constant function of input voltage V x . Illustrations of the 
quantization function are shown in Fig. B-l for two-level and three-level quantization, with 
quantized output voltage V x plotted as a function of input voltage. Discriminator levels that 
determine quantization boundaries are placed at the voltage values that cause the minimum SNR 
loss, based on the assumptions of low single-sample SNR (e.g., <0.1) and Gaussian noise. For 
two-level sampling, the single discriminator level is at zero voltage. For three-level sampling, 
the two required levels are approximately -0.61 G^ and +0.61 G^, as illustrated in Fig. B-l and 
discussed in Subsection B.l. 3, where G^ is the variance of either noise component. 

In the presence of Gaussian noise with a standard deviation of G^ for each quadrature 
component, the expectation value for the quantized voltage at a particular sample point can be 
computed as 


<V X > T1 = Q(V x )P(V x ,S x )dV x (B.3) 


where < > T) denotes the ensemble average over noise, and where the probability density function 
is given by 




1 Q-(V x -S x ) 2 /2al 

i2n Gt) 


(B.4) 
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Figure B-l. Quantization Function for Two-level and Three-level Quantization 
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in which S x is the value of the real component of the signal at the particular sample point under 
consideration. For the low single-sample SNRs, Sx/c^ , expected for the TR receiver, the 
expectation value of a quantized sample is linear with respect to the signal S x and can be 
approximated by 




(B.5) 


where the proportionality factor q M is specific to the quantization function. Values for q M are 
derived below for two-level and three-level quantization. 

The variance of the noise on the quantized voltage is defined by 

^.<(V,-<V X >,) 2 >, (B.6) 

which, in the case of low SNR, can be approximated by 

a ! ~ a iT = <^lx >n (B.7) 

where r| x is sampled noise. Thus, the variance of the noise on the sampled voltage is 
approximately equal to the variance of the noise sampled without signal and is given by 


a~ = Q 2 (V x )P(Vx,0)dV x (B.8) 


where S x has been set to zero. 


The single-sample voltage SNR after quantization is defined by 


SNR s = 




which is approximately given by 


SNR s - 




(B.9) 


(B.10) 


in the low-SNR limit. Expressions for SNR are developed below for two-level and three-level 
quantization. 


B.1.2 Two-Level Quantization 

The sampler converts the filtered baseband signal to a two-level (1-bit) form, with a 
quantization function defined by 
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Q(V) = +1, V>0 
Q(V)= -1, V < 0 


(B.ll) 

(B.12) 


This quantization algorithm is sometimes referred to as sign sampling since the single bit 
representing the sample is the sign of the voltage at the sample point. When system noise greatly 
exceeds the signal, two-level sampling causes negligible distortion on average, of the signal 
component of the sampled voltage, but does decrease SNR by approximately 20%, as outlined in 
the next paragraph. 

For two-level quantization, Eq. (B.3) reduces to 


< V x >q — 



(B.13) 


If the SNR is low (i.e., Sx/oy, < 1), this integral can be expanded in the form 


< Vx >n - 



(B.14) 


When SNR is so low that the S 3 term can be ignored, the sampled signal will be given on average 
by the linear form of Eq. (B.5) with 



(B.15) 


Since the quantization level is either +1 or -1, the variance of the sampled noise computed in 
Eq. (B.8) becomes 



(B.16) 


Substitution of these results in Eq. (B.10) shows that the single-sample SNR after sampling is 
related to the SNR before sampling, S x /CTn , by 


SNR s - ^ 

V it 


(B.17) 


Thus, the voltage SNR is decreased by approximately 20% by two-level sampling. 

In the strongest P channel, Ll-P, the largest analog SNR, S x /oV|, is approximately 0.17 
near maximum antenna gain (about 7 dB with the choke-ring antenna). Thus, the cubic term in 
Eq. (B.14) contributes less than 0.5 % to the sampled signal. Tests and theory indicate that this 
worst-case cubic term leads to negligible errors in phase and to less than 1-cm errors in P delays 
If necessary, the S 3 delay error could be further reduced on the basis of concurrently observed 
SNR. 


For the C/A channel, S x AXq can be as large as 0.25 for the TR receiver, so that the cubic 
term becomes relatively large (e.g., 1% of the linear term). A calibration attempt has suggested 
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that the worst-case S 3 error in delay for C/A delay is less than 2 cm, but the effect was too small 
to measure in the attempted scheme. Again, the S 3 error in phase is negligible. In principle, a 
correction based on concurrently observed SNR could be applied to reduce the S 3 delay error, but 
such a correction is not applied in the current software. 


B.1.3 Three-Level Quantization 

When an analog signal is subjected to three-level sampling, the quantization function is 
given by 

Q(V) = +1, V > D (B.18) 

Q(V)= 0, |V| <D (B.19) 

Q(V) = -1, V < -D (B.20) 

where D is a positive discriminator parameter. For this quantization function, the expectation of 
the sampled voltage in Eq. (B.3) reduces to 


<V X >t| 



(B.21) 


where 


Zu = 


D-Sx 


(B.22) 


and 


z L = 


-D-Sx 


(B.23) 


Calculation of the first, second, and third derivatives of Eq. (B.21) with respect to S/an 
leads to the expansion 

<V x >n = e -” 2/2 

where 

u s (B.25) 

an 


Sx (l-u 2 )/S x \ 3 


(B.24) 


Thus, the proportionality factor in Eq. (B.5) becomes 



(B.26) 
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Note that Eq. (B.24) reduces to Eq. (B.4) for u = 0, as one would expect. For the discriminator 
level that corresponds to u = 0.612, Eq. (B.24) becomes 


<%> n * 0.662 



(B.27) 


which indicates that the proportionality factor becomes q M = 0.662. 

Comparison of Eq. (B.27) with Eq. (B.14) indicates that the linear term is 17% smaller 
for three-level sampling than for two-level sampling. As shown below, this loss in amplitude is 
more than offset by a decrease in RMS noise. In addition, the cubic term is a factor of 1.6 less 
corrupting for three-level sampling, which indicates that the resulting phase and delay errors will 
be correspondingly lower. 

For three-level sampling, the variance of the sampled noise in Eq. (B.8) reduces to 



—2— f e’ z2/2 dz 

V2¥ I 

Ju 


(B.28) 


When u = 0.612, one obtains 


o~ = 0.5405 

n 


(B.29) 


By combining the linear term in Eq. (B.27) and this estimate of noise, as prescribed in 
Eq. (B. 10), one obtains the single-sample SNR 


SNR s = 0.900 (B.30) 

<Jn 

This result shows that three-level sampling of an analog signal causes a voltage SNR loss of 
approximately 10% in the case of low SNR. Thus, as indicated by Eqs. (B.17) and (B.30), three- 
level sampling provides about 13% better voltage SNR (approximately 1 dB) than two-level 
sampling. 

By using the above expressions, one can show that the assumed discriminator level of 
u = 0.612 minimizes SNR loss. Form the single-sample SNR by extracting the leading term of 
Eq. (B.24) and dividing by the square root of Eq. (B.28). The resulting function is maximized 
when u is approximately equal to 0.612. 


B. 1.4 Highly Digital Front End 

In the highly digital front end presented in Appendix A, two-level sampling is carried out 
at RF at a very high rate (about 450 MS/s for LI and about 490 for L2). The samples are then 
subjected to a sum-and-dump operation to reduce the rate to 20.456 MS/s and increase SNR. 

The sum, which covers 22 points for the LI channel and 24 for L2, increases the number of bits 
required to represent a sample to a relatively high value (6 bits). To keep a minimum-bit design 
in the digital baseband processor, the sum output values can be requantized to three levels, -1,0, 
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and +1, with the integer discriminator level set as closely as possible to the value that minimizes 
SNR loss. As mentioned above, the optimum discriminator level settings are approximately 
-0.61 a and +0.61a for ordinary three-level sampling of an analog voltage. Integer discriminator 
levels can be set in analogy with ordinary sampling. Given the RF-filter bandwidths in Fig. A-l, 
for example, noise is uncorrelated between samples, the RMS value for noise leaving the sum is 
VN , where N is the number of points (22 or 24) in the sum. As an approximation, therefore, the 
requantization discriminator parameter can be computed as 0.612 VN, which equals 2.871 for LI 
and 2.998 for L2. These values are both rounded to 3, the nearest integer, which leads to 
discriminator values of -3 and +3 for both the LI and the L2 channels, given the implementation 
in Fig. A-l. Since the number of points in the sum is even for both LI and L2, all sum values are 
even, and only the sum values of -2, 0, +2 are requantized to zero; all other sum values become 
either -1 or +1. 

The cubic effect introduced by the initial two-level sampling at RF is very small for both 
channels due to the very low single-sample SNR caused by the very wide noise bandwidth (e.g., 
450 MHz) entering the RF sampler. Thus, even though the initial two-level sampling degrades 
the SNR by 20%, the initial 1-bit samples are, on average, essentially linear with respect to the 
input signal. To first approximation, the N-point sum, on average, preserves this linearity. As an 
approximation, therefore, the effects of integer requantization can be modeled as a three-level 
quantization of an analog signal with an SNR equal to the SNR of the integer samples leaving the 
sum. This model should closely approximate the amplitude and noise effects of requantization. 

When the effects of three-level requantization are combined with the effects of the initial 
two-level sampling at RF, therefore, one obtains the following approximate expression for the 
SNR for each 20.456-MHz sample: 


SNR s ~ 0.900 (B.31) 

V 71 

where S x /<7r| is the single-sample SNR that would have been obtained by sampling at 20.456 

MHz with an "infinite" number of quantization levels, V2/7t is the loss due to the initial two-level 
sampling, and 0.900 is the loss caused by three-level requantization. The combined effect 
becomes 


SNR s ~ 0.72— (B.32) 

Gt| 

Thus, the combined SNR loss due to both quantization steps in the example digital front end is 
about 28% in voltage SNR or 2.9 dB. If the bandwidth of the filter preceding the RF sampler is 
greatly reduced (e.g., from 450 MHz to 90 Mhz) so that the signal is greatly oversampled, much 
of the SNR loss caused by the initial two- level sampling can be "regained." 

In analogy with the initial two-level quantization, a cubic term is introduced by the three - 
level requantization. An approximate analysis will be carried out through analogy with the 
results presented in the preceding subsections. Since the single-sample SNR at requantization is 
about the same as the single-sample SNR at the point of two-level quantization in the analog 
front end (see Subsection B.1.3), the S 3 effect is about 1.6 times smaller for the digital front end 
due to the smaller coefficient of the cubic term at requantization (compare Eqs. B. 14 and B.27). 
Thus, the S 3 error is estimated to be about a factor of 1.6 smaller for the digital front end than for 
the analog front end (see Subsection B.1.2). 
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B.2 Three-Level Quantization of Counter-Rotation Sinusoids 

The model phasor of Eq. (4.2) is quantized on each component to three levels in the 
TurboRogue receiver. As illustrated in Fig. B-2, the quantization consists of amplitude values 
equal to -1, 0, or +1 with subcycle divisions of 1/8, 3/8, 1/8, 3/8 cycle. This quantized form for a 
sinusoid greatly simplifies the hardware but causes a slight increase in SNR. If one Fourier 
decomposes the three-level sinusoid, the coefficient of the fundamental (i.e., the desired 
counterrotation sinusoid) is 1.176, which is 17.6% higher than for an unquantized sinusoid. This 
increase in amplitude is canceled by greater effective noise, which increases by a factor of VI. 5 
(see Subappendix C.2.1). The net effect is a decrease in voltage SNR of about 4% (-0.35 dB) 
due to three-level quantization. Greatly reduced complexity in the hardware justifies this small 
sacrifice in SNR. 

Higher harmonics are introduced by three-level quantization of the counterrotation 
sinusoids. If the baseband carrier frequency is sufficiently high, these unwanted harmonics are 
reduced to a negligible level by the filtering effect of the correlation sum. A high baseband 
frequency is produced by offsetting the downconversion frequency. For example, the LI 
baseband frequency is on the order of 300 kHz, and its third harmonic is on the order of 900 kHz. 
This harmonic is counterrotated to about 600 kHz by the baseband processor. A 20-ms 
correlation sum would then reduce its amplitude by a factor greater than 2*600 kHz * 0.02 s = 
24,000 relative to the primary signal. At this reduced level, the third-harmonic term would 
introduce a phase error less than 7 [icycles. Even greater reduction in error will result from 
averaging over longer time intervals (e.g., 1 s). Thus, harmonics introduced by three-level 
quantization are reduced to negligible levels in measured phase. It is assumed that the incoming 
signal has no tones at the harmonic frequencies. 
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APPENDIX C 


NOISE-CORRELATION ANALYSIS 


C.1 Correlation of Noise Between Samples 

Let Rn(x) be the analog autocorrelation function for each of the two components of the 
complex noise, T|, on the baseband signal in Eq. (3.2). If the receiver in the ideal case is 
characterized by a rectangular bandpass with a two-sided bandwidth equal to W, the analog 
autocorrelation function of the noise produced by double-sideband downconversion will be given 

by 


R T1 (t) = 


sinpcWx] 

JtWx 


(C.1) 


after normalization to one at x = 0. For this ideal configuration, the sample rate f s is equal to W 
in a Nyquist implementation, and the zeros of the analog autocorrelation fall at sample-interval 
separations of 8x = 1/W = l/f s . Thus, the ideal autocorrelation function is equal to one at the 
origin and zero at other sample spacings, indicating that noise on a given sample is uncorrelated 
with the noise on any other sample. The actual bandpass is not rectangular, but the system 
implementation is sufficiently close to ideal to allow this model to serve as an approximation for 
noise analysis. 

The autocorrelation function for two-level sampled noise (van Vleck and Middleton, 
1966) can be expressed in terms of the analog autocorrelation function as 

Rti ft) = f sin- 1 [Rn (x)] (C.2) 


This relation shows that the digital autocorrelation is equal to one at the origin, with nulls at the 
same points as in the analog case. Thus, the noise is still uncorrelated between different samples 
after two-level sampling. A similar conclusion applies to three-level quantization. 

As analyzed in Appendix D of this document, the two quadrature noise components are 
assumed to be uncorrelated with one another and to have the same noise characteristics. With 
these assumptions and the features of the digital autocorrelation functions described above, one 
can easily show that the sampled complex noise satisfies the relations 

< il(tj ) r\ *(tk»T, = 2o~ 5jk (C.3) 


<Tl(tj)tl(t k )> 11 = 0 (C , 

where the tilde denotes sampling quantization of each component and G~ is the RMS value of 
sampled noise for each component. Expressions for G~ can be found in Eqs. (B.16) and (B.29) 
for two- and three-level quantization, respectively. 
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C.2 Code Mode 


C.2.1 Two Complementary Correlations for the Noise on the Correlation Sums 

Two complementary correlations are derived to begin assessment of the noise on the 
complex correlation sums. The first calculation determines the complex interlag correlation 

Pl(h O = <(U/ - <U/>)(Uj*-- <U;'>)> 
and the second the complementary correlation 

p 2 (/,/’) = <(U/ - <U; >)(U/’- <U/'>)> 

where < > indicates an ensemble average over both noise and code. 

For weak signals, the first correlation is approximated by 

Pl(U')«<u/ U;.> 

Substituting the definition of correlation sum in Eq. (4.5), one obtains 

pi(/, n « X <v(tj)v(t k ) v m (tj,/ )v^(t k , /•) > 

j.k 

where, for brevity, t k refers to t R + ks. For a small single-sample SNR, the sampled voltage is 
mostly noise, and Eq. (C.8) becomes 

P 1 (A l ) 5=5 X ^ "nCtj ) R (tk) >n ^ V m (tj , l ) V m (t k , / ) >p (C.9) 

j.k 


(C.5) 

(C.6) 

(C.7) 

(C.8) 


where the noise and code averages have now decoupled. In the case of ideal sampling, Eq. (C.3) 
reduces this equation to 

Pl(j>0~ 2 g~X <V m (t k ,/)Vm(t k ,/ , )>p (C.io; 

k 

When the correlator model in Eq. (4.1) is substituted, one obtains 


~ 2g~X |C( ( Pm)| 2 <P(t k -'t m -/ s) P(t k -t m -/'s)>p 

k 


(C.11) 


The expectation value in this equation can be expressed in terms of the P-code 
autocorrelation function as defined in Eq. (2.12), which yields 

Pi (hi ) ~ 1 Q CPm) 1 Rp(t k -x m -/s, ls-l s) 

k 


(C.12) 
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For three-level quantization, one can easily show that the magnitude of the quantized counter- 
rotation phasor is equal to 


|Q(Pm)| 2 =1.5 + square wave (C. 1 3) 

where the square wave has a frequency equal to four times the carrier rate. In the sum over time, 
this square wave with a frequency on the order of 1.2 MHz is typically orthogonal to the 
autocorrelation function with its frequency of about 10.23 MHz, particularly when Doppler 
shifting is incorporated. Thus, Eq. (C.12) becomes 

Pj(/, /’ ) = 3a|X R p( t k-'Cm-/s, /s-/*s) (C.14) 

k 

As indicated in Eq. (2.14), the time average of the P-code autocorrelation function is a triangle 
function, so that Eq. (C.14) becomes 

Pt(^> l' )~ 3N s ct~ A(Z - /' ) (C.15) 


where N s is the number of samples in the correlation sum, and A is the triangle function shown in 
Fig. 2-3. 


If the counterrotation phasor were not quantized, the phasor magnitude in Eq. (C.13) 
would be equal to one for all time points and the coefficient in Eq. (C.15) would be 2 rather than 
3. This comparison shows that three-level quantization of the counterrotation sinusoids increases 
the variance of the sum noise by a factor of 1 .5. 

Analogous to the preceding analysis, the second interlag correlation in Eq. (C.6) reduces 
to 


P2(U' )*2 < T l(tj) r l( t k)>T|< V m(tj,/ )V m (t k ,/' )>p 

i.k 


(C.16) 


Because of Eq. (C.4), one obtains 


P2(U’) = 0 


(C. 17) 


which shows that this correlation is zero for all lag combinations. 


C.2.2 Interlag Correlations and RMS Sum Noise 

Let x; and y / be the real and imaginary components of the complex sum noise, u / - <u/>. 
Eqs. (C.5) and (C.15) show that averages of the products of these components satisfy the 
following two relations: 


< x/ x/>+ <y/ y /•>« 3N S a~A(/ ) 


(C.18) 


<x/-y /> - < x/ yr> = 0 


(C.19) 
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while Eqs. (C.6) and (C.17) show 


<x z x/>- <yz yr>~ 0 

(C.20) 

< x/ ■ y/ > + < x/ yz > ~ 0 

(C.21) 

Eqs. (C. 18) and (C.20) can be combined to yield 


< xz xz ■> = < yz yz ■> ~ 1.5 N s g~A(/ ) 

(C.22) 


while Eqs. (C.19) and (C.21) show 


<x/ y;> = 0 


(C.23) 


This last equation indicates that the noise on the complex correlation sum is uncorrelated 
between real and imaginary components, regardless of lag. 


For l = l ', Eq. (C.22) leads to the result that the variance of the noise on either component 
of the correlation sum, which is defined by 

a^= <x ; 2 >= <y z 2 > (C.24) 


is given by 

a 2 =1.5 a~N s 


(C.25) 


and is therefore the same for all lags. The RMS sampled noise is given by either Eq. (B.16) 
or Eq. (B.29), depending on the number of quantization levels. N s is the total number of samples 
in the correlation sum. As mentioned above, the factor of 1 .5 represents the increase in noise 
variance due to three-level quantization of the counterrotation sinusoids. 


For either component, the normalized correlation of noise between lags is obtained by 
dividing Eq. (C.22) by Eq. (C.25), which gives 


J-<x/ x/'>--L<y, y r > = A (/-/') 

as as 


(C.26) 


Thus, the normalized correlation of the sum noise between lags is given by the triangle function 
obtained by autocorrelating the pure P-code (see Subsection 2.3). Consequently, the normalized 
correlation of noise between adjacent lags is approximately 0.5; for a lag separation of 2 or 
greater, the correlation is nearly zero. 
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C.3 P-Codeless Mode 

C.3.1 Two Complementary Correlations for the Noise on the Correlation Sums 

The two complementary correlations in Eqs. (C.5) and (C.6) derived above for code- 
mode operation must be rederived for P-codeless operation. Based on Eq. (5.1), the correlation 
in Eq. (C.7) becomes 

pi(U') - <X V,(tj -/s) V 2 (tj) Qtj) Vi*(t k -/'s) V 2 (t k )C(tk)> (C.27) 

j.k 

where the ensemble average is over noise on LI and L2 and over the P code. For brevity, t k 
refers to tR + k s and l represents L + /. The counterrotation phasor, C, has now been explicitly 
labeled with the appropriate time argument instead of the corresponding phase argument. Since 
noise dominates in the limit of low single-sample SNR, this expression becomes 

pi (/,/’) «<XBi(tj ^s)nfo)Cty)Th( t k-/ , s) Tte(t k ) C'ctk)^ (C.28) 

j.k 


where the ensemble average is now only over noise. Separating the ensemble averages for LI 
noise and L2 noise leads to 

PiC/,/’) -S <ni(tj-/s)^r(tk-/ , s) > T | 1 < : H 2 (tj)B 2 (tk)> ll2 C(t j )C*(t k ) (C.29) 

j-k 

By substituting Eq. (C.3) twice with appropriate subscripts, one obtains 

pi(/, /’) » X 2 <4 5 K k -r 2c j| 2 8k, j C(tj) C*(tk) (C.30) 

j>k 

Because of the Kronecker delta functions, this expression is nonzero only when l is equal to l \ 
and in that case, the double sum collapses to a single sum: 

Pl(U’) ~48 M .tf^ 2 X|C(t k )| 2 (C.31) 

k 

If one uses Eq. (C. 13) and assumes the square wave averages to zero, one obtains 

pi(/,Z’) ~ 6N s a~j a~ 2 S/,r (C.32) 

where N s is again the number of points in the correlation sum. 

In analogy with the derivation in Subsection C.2.1, one can easily show that the second 
correlation, p 2 , defined in Eq. (C.6), is equal to zero for all lag combinations, as in Eq. (C.17). 
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C.3.2 Interlag Correlations and RMS Sum Noise 

By using Eqs. (C.17) and (C.32), and following the derivation in Subsection C.2.2, one 
can show that the components of the complex sum noise for L1*L2 satisfy the following 
relations: 


<x/ x/> = <y/ y/> * 3 Nsa-j a~ 2 5/, r 


(C.33) 


and 


<x/ y/'> = 0 


(C.34) 


This last equation indicates that the noise on the complex correlation sum is uncorrelated 
between real and imaginary components, regardless of lag. Equation (C.33) shows that the noise 
on the real and imaginary components of the correlation sums is respectively uncorrelated 
between lags and that the variance of the noise on either component of a correlation sum is given 
by 


gI 


3NcC7 


m 



(C.35) 


Again, the RMS sampled noise <j~ is given by either Eq. (B.16) or Eq. (B.29), depending on the 
number of quantization levels. As in Eq. (C.25), this variance is increased by a factor of 1.5 by 
the three-level quantization of the counterrotation sinusoids. A factor of 2 is introduced by 
correlating one quadrature signal with another, as opposed to correlating a quadrature signal with 
a model code. 

For either component, the normalized correlation of noise between lags is obtained by 
dividing Eq. (C.33) by Eq. (C.35), which gives 


J-<x z x/>« ^-<y/ y /•>« 5 /,/■ 
Gu Gu 


(C.36) 


These results allow a comparison of sum-noise correlations in the P-codeless mode and 
code mode. Equation (C.36) shows that sum noise is uncorrelated between any two separate lags 
in the P-codeless mode. In the code mode, on the other hand, Eq. (C.26) shows that sum noise is 
uncorrelated between any two lags with a separation greater than one sample but is strongly 
correlated (0.5) between adjacent lags. 


C-6 






APPENDIX D 


QUADRATURE NOISE STATISTICS 

This appendix analyzes the statistics of the noise on the two quadrature components of 
the baseband signal, including the correlation between those components. 

Analogous to Eq. (2.17), the noise at RF before filtering can be written as 


N R (t) = 



A 11 (v)e 2lcivt dv + c.c. 


(D.l) 


where c.c. denotes the complex conjugate of the first term, v is the RF frequency (i.e., fL + f in 
Section 2), and is the Fourier distribution of the noise at RF before any receiver filtering. 
Note that v instead of f is now used as the independent frequency variable. In the front-end 
processing that takes the signal to baseband, the noise component is subjected to the same 
filtering and do wncon version as the signal, as modeled in Subsection 3.1. After front-end 
processing, the (real) RF noise in Eq. (D.l) is transformed to complex baseband noise given by 


q(t) = 



G r (v) At,(v) e 2K '[vt- cp h ] dv 


(D.2) 


where 9h is downconversion phase and Gr is the complex system filter for the receiver. The 
same operations that transformed Eq. (2.19) to Eq. (3.2) have been applied to the noise in 
Eq. (D.l), except for the antenna filter. This equation implicitly assumes that sine 
downconversion phase is exactly 90° out of phase with the cosine phase. An offset error in sine 
phase can be modeled as described below. 

The cosine component of the baseband noise in Eq. (D.2) can be modeled as 


if 

Jo 


T| c (t)= 4-1 G R (v)A 11 (v)e 2,ti ( vt - < Ph) dv + c.c. 


Similarly, the sine component of the baseband noise can be modeled as 


(D.3) 


■i 

Jo 


r| s (t) = \ I G R (v)A T1 (v)e 27Ci ( vt - < Ph+ «) dv + c.c. 


(D.4) 


where 0 is a phase offset ideally equal to 0.25 cycles when the downconversion phase for the 
sine component is exactly 90° away from the downconversion phase for the cosine component. 

The normalized cross-correlation function for the two components is defined by 


R x (t-t’ ) = -V < "He (t) Tls(t' )>^ 

0V| 


(D.5) 
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where it is assumed that the noise power is the same for the two components and where is the 
RMS noise for either component. (In practice, the RMS noise can be different for the two 
components, but such a difference has no effect on the final result from code correlation as long 
as the receiver amplitude-versus-frequency and phase- versus-frequency profiles along the two 
component paths are the same. This insensitivity to relative gain is a consequence of two-level 
sampling, which makes the per-sample SNR the important quantity in TR digital signal 
processing.) If Eqs. (D.3) and (D.4) are substituted in Eq. (D.5) and the averaging operation and 
the integrals are interchanged, one obtains 


R x (t-t') = 



Gr(v)G r (v')<A 11 (v)A t1 (v' )> t) e 27ti (vt-<ph+v't'-<ph+ 0) dv dv' 



Gr(v)Gr(v')<A t1 (v)A* (v' )> n e 27ti (vt-<Ph-v't’+ <ph-0) dvdv' 


(D.6) 


where c.c. denotes complex conjugation of both of the first two terms. If the RF noise is 
stationary, one can show 


<A 11 (v) A^v' ^ = 0 when v, v' >0 
and 

<A T1 (v)A*(v’)> T1 =i-N 0 (v)8(v-v’) 

where No is the noise power spectrum at RF before receiver filtering. With this definition of No, 
the total noise power at RF, including both positive and negative frequencies, is NoW, where W 
is the full noise bandwidth. 

When Eqs. (3.1), (D.7), and (D.8) are substituted in Eq. (D.6), one obtains 


Rx(t-t')= -t- | G R ( v ) | 2 N 0 (v ) cos[2 7t (v - f h ) (t - 1 1 ) -0] dv (D.9) 

4 <*n Jo 

The use of the linear form of downconversion phase in Eq. (D.9) assumes that, over a correlation 
interval, nonlinear time variations in downconversion phase are negligible, so that the 
"drift," 9hd, disappears in the difference even when t is not equal to t’. Define the quadrature 
phase error A0 as the deviation of the quadrature phase shift 0 from 90°: 

0 = 0.25 + A0 (D.10) 

Eq. (D.9) can now be rewritten in terms of this phase error as 


(D.7) 


(D.8) 
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(D.ll) 


Rx(t-f) = -L- \ | G R ( v ) | 2 N 0 (v ) sin[2 K (v - fh) (t - 1') +A9] dv 

4 tfJ o 

When the quadrature channels are exactly 90° out of phase ( A9 = 0), Eq. (D.l 1) shows that the 
two quadrature components will be uncorrelated for all values of t and t’ if and only if the noise 
spectrum after receiver filtering is symmetric about the downconversion frequency. For 
example, if the noise spectrum is flat and the composite filter of the receiver is symmetric about 
fh, the two quadrature components will have independent statistics, and processing both 
components will improve SNR by 3 dB relative to processing only one. 

A similar analysis shows that the variance for either of the baseband quadrature noise 
components is given by 


<y2 s <q c 2 > = <ni > = 1 f |Gr(v)| 2 N 0 (v) dv (D.12) 


When the spectrum after receiver filtering can be approximated by a rectangle with a double- 
sided noise bandwidth of W, the noise variance on each quadrature component becomes 

^ = 1|Gr| 2 N„W (D.13) 

The factor of 1/4 (or 1/2 for a^) arises from quadrature-component extraction in Eqs. (D.3) and 
(D.4) and is of no ultimate consequence since the same factor would enter both the signal and the 
noise. 


Given Eq. (D.13) and the assumption of a rectangular spectrum, the correlation between 
quadrature components in Eq. (D.l 1) becomes 


R x (t-t') = 


sin[JcW(t-t')] 

TtW(t-t') 


sin(A9) 


(D.14) 


where it is assumed that the downconversion frequency fh is at passband center. This result can 
be used to place an upper limit on quadrature phase error. For the worst case of t - 1' = 0, the 
correlation is equal to A9 (in radians) when A9 is small. For example, to keep intercomponent 
correlation below 0.05, the quadrature phase error must be less than 0.05 radians (about 3 
degrees). 

When the digital front end discussed in Appendix A is analyzed, one finds by a similar 
analysis that the correlation between the noise on the quadrature components can be reduced to a 
negligible level by making the composite bandpass, including digital filtering, symmetric about 
the effective downconversion frequency. However, because the quadrature components are 
generated through sampling rather than by downconversion through use of an LO signal, 
quadrature phase error, A9, is not a problem in the digital front end. 
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APPENDIX E 


THE C/A CORRELATION FUNCTION 


Derivation of the C/A correlation function closely parallels the derivation for P. 
Differences are a result of the cyclic nature of the C/A code. The spacecraft generates the C/A 
code at a rate of 1.023 x 10 6 c-chips/s by filling every millisecond with the same sequence of 
1023 chips, where each 1-ms sequence begins at an integer multiple of 1 ms, GPS satellite time 
The ideal unfiltered C/A code can be modeled as 


C(t)= I C s (t-t-)m(t')dt' (E.l) 


where C s (t) is the sequence of 1023 chips. C s (t) can be represented by a form similar to 
Eq. (2.1): 

512 

Q(t)= I bjB’Ct-tj) (E.2) 

j = — 512 

where B 1 is a function like that in Fig. 2-2, but with a width of T c = (1.023) 4 (is; where the 
coefficients bj are the ostensibly random chip values of -1 or +1, and where 

tj = j T C (E.3) 

Note that the center chip, 513, has been placed at tj = 0. The sampling function is defined by 


4- oo 


m(t)= £ 8(t-(k + 0.5)Ti) 


k=- 


(E.4) 


where T i = 1 ms, and 8(t) is the Dirac delta function. 

Paralleling the P-code derivation, the Fourier transform of the ideal C/A code over an 
accumulation interval of width T and center t is given by 


A(f) 


/•t+T/2 

"l 

/ t-T/2 


C(t)e-2rci ft dt 


(E.5) 


The Fourier transform of the single 1023 sequence is given by 

fTi/2 

A s (f)= C s (t)e-2*iftdt 

J- Ti/2 


(E.6) 


which becomes 
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A s (f) = T c 5!S^£bje-2* if >i 


(E.7) 


where tj is at the center of the jth chip given by Eq. (E.2). The Fourier transform of the sampling 
function is given by 


{ 


rt+T/2 

W(f) = | III(t)e~ 27tift dt 
t-T/2 


(E.8) 


which becomes 


W(f ) = e 


„2wift sin(NiTtfTi) 


(E.9) 


sin(7tfTi) 

where Ni = T/Tj is the number of 1023 sequences repeated over the transform interval T. 

Based on the convolution theorem, the Fourier transform of the C/A code is given by 

A(f) = A s (f)W(f) (RIO) 

so that the power spectrum becomes 


|A(f)| 2 = |A s (f)p 


sin(Ni7tfTi) 


sin(7tfTi) 


(E. 11) 


In the limit of large Ni, one can show 


sin(Ni7tfTi) 


sin(7tfTi) 
so that the power spectrum becomes 


-^2 8 ( f - n/T i) 


(E.12) 


|A(f)| 2 «|A s (f)| 2 ^-X 5 ( f - n/T i) 

T 1 n 


(E.13) 


where A s (f) is given by Eq. (E.7). Thus, the C/A power spectrum is approximately equal to the 
power spectrum of the 1023 sequence multiplied by a sum of delta functions that select points at 
multiples of 1 kHz. 

The ideal, unfiltered C/A signal at RF can be represented by 

Si(t) = C(t) sin(27tf L t) (E.14) 

From this point, modeling of spacecraft and receiver-front-end processing parallels the analysis 
of Subsection 2.4 and Section 3 and will not be repeated here. 
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Correlation analysis, however, should take into account "code noise", defined as the 
deviation of results from the C/A code from the results from a long random sequence. Unlike the 
P code, the short repetition interval and nature of the C/A code makes the long random-sequence 
approximation fairly inaccurate for the C/A code, at least with respect to the amplitude-vs-lag 
dependence of the correlation function. Since the C/A code has only 1023 separate chips for any 
correlation interval greater than 1 ms, one would expect the code noise effect to be much larger 
for C/A than for P. For example, if the C/A chips were independent of one another, C/A code 
nois e woul d theoretically cause a fractional amplitude error at the early or late lag on the order of 
0.5 /V 1023 = 0.015 , on average, relative to peak amplitude, given an accumulation interval 
greater than 1 ms and 500-ns lags. (In practice, since the C/A code is not a truly random 
sequence, the magnitude of this correlation turns out to be approximately equal to either 0.03 or 
0.0007, depending on the C/A sequence.) On the other hand, fo r the P co de with 50-ns lags, the 
corresponding fractional amplitude error is on the order of 0.5 /V 200,000 = 0.001 , given an 18- 
ms accumulation interval and 50-ns lags. Since the maximum voltage SNR obtained by TR for 
an 18-ms interval is about 100, the minimum fractional amplitude error due to system noise is 
about 0.01 . Thus, while P-code self noise is very small compared to system noise and is 
therefore negligible, C/A-code self noise potentially can be an important error source. 

To obtain a correlation function for C/A that accounts for code noise, one can perform an 
analysis similar to the analysis in Section 4, including an ensemble average over noise, but 
omitting the ensemble average over code. Such an analysis yields the following approximation 
for the C/A correlation function: 


( + °o 

G(f) 


LMQi! e - 2*i [f(t, - lb) - <p.] df 

TT C 


(E.15) 


where b is the lag-offset interval for C/A and 


A^ N ,T c c,3* Sh,M 

°n 7l<PdT 


(E.16) 


and where the power spectrum | A(f)| 2 is given by Eq. (E. 11). As suggested by Eq. (2. 11), the 
product, TT C , dividing | A(f)| 2 in Eq. (E.15), normalizes the power spectrum so that it corresponds 
to a (sinx/x) 2 function that equals 1.0 at the origin. Note that this equation is the same as 
Eq. (4.22) for the P code if the normalized power spectrum corresponds to a (sinx/x) 2 function. 

Substitution of Eq. (E.13) in Eq. (E.15) yields 


u/(t R )«^e 2 *W£G(f„) 


l A s(fn)l 2 e - 2 Ki [f„ Ct r - lb) - cp,(f„)] 
TlT c 


(E.17) 


where f n = n/Tj and where the satellite-specific "1023" Fourier component A s is given by 
Eq. (E.7) in terms of specific C chips. The power spectrum is again normalized to correspond to 
a (sinx/x) 2 function that equals 1.0 at the origin. Thus, the C/A correlation function can be 
computed as a sum over frequency points with 1-kHz spacing. The factor of 1/Ti,which is equal 
to the 1-kHz separation between frequency points, can be viewed as a replacement for df in the 
conversion from integral to sum. 
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Based on the preceding equations, one can use each satellite-specific 1023 code sequence 
to accurately model the effects of C/A code noise. For example, if one uses Eq. (E.13) to obtain 
the power spectrum of satellite 1 based on 1-ms of code, one obtains the profusion of noisy 1- 
kHz points shown in Fig. E-l. Based on the normalization for |A s (f)| 2 shown in Eq. (E.17), this 
spectrum has been normalized so that the power would approach 1 .0 at zero frequency, were it 
not for code noise. A more informative picture is gained by averaging this noisy spectrum over 
larger frequency bins. For an averaging interval of 35 kHz, for example, one obtains the spectra 
plotted in Figs. E-2(a), E-2(b), and E-2(c) for satellites 1, 15, and 22, respectively. These 
satellite examples each correspond to one of only three values for the code-induced fractional 
amplitude errors found at the early and late lags, namely, the Butterworth-filtered values 
approximately equal to -0.0007 (satellite PRNs 1, 2, 3, 4, 5, 6, 9, 10, 1 1, 12, 13, 14, 16, 18, 19, 
20, 23, 25 through 32); +0.0315 (7, 15, 17, 21, 24); and -0.0318 (8, 22), with all based on "500- 
ns" lags. (Alternately expressed, the correlation of adjacent chips across 1023 chips is 
approximately equal to twice the respective values just quoted.) 

The spectra in Figs. E-2(a), (b), and (c) can be substituted in Eq. (E.17) to produce the 
amplitude-versus-delay responses plotted in Figs. E-3(a), (b), and (c). To better display the code- 
induced amplitude error, however, these responses are each plotted relative to the code-noise-free 
response (i.e., long-random-sequence approximation) in Fig. 4-3(a). The three aforementioned 
fractional amplitude errors are evident at delay values of ± 0.25 c-chip ("500 ns"). 

For TR, C/A code noise is not a significant error source for phase and delay 
measurements. Code noise has a negligible effect on phase since phase is extracted from the 
prompt lag, where signal and model are very nearly aligned to eliminate code sign. Furthermore, 
given the early-minus-late algorithm for extracting C/A delay, as shown in Eq. (6.7), code noise 
has no effect on delay either, since it changes the amplitude of the early and late lags by the same 
amount. Thus, although code noise is potentially an important error source for C/A, it causes no 
degradation in phase and delay when TR uses the processing algorithms presented in this report. 


E-4 



NORMALIZED SPECTRAL POWER 


6 


5 


4 


3 


2 


1 


0 


r - — i 1 1 

o 

o 

1 1 1 1 

1 1 1 1 

1 I I l 

— o 


o = satellite 1, 1-kHz S] 
line = (sinx/x) 1 2 

racing ~ 

J- 

o 

0 o 

“ o 

” o o 0 



- 

_ <p 

St. o o ° 0 

e 

XlO 0 Q 

* o 8 Ors ° 

OO °O rt 

o o 0 ° O O 



- 

- 7-'7“o* — 

3 ° ° °* ° .. 

- ®b ®o o° ° o c 

’ Q°°a ° ° °%>0 OO ® o ° o 

“ 8 o ft ° O® o° ° 

rv n an O Qnff ^ 

> <? 
o 

OO 


- 

|MHN 

•••o’ 

O 

1 

O 



0 


0.5 


1 1.5 

FREQUENCY, MHz 


2 


Figure E-l. Power Spectrum of the C/A Code for Satellite 1, with 1-kHz Spacing 





Figure E-2. Power Spectra of the C/A Code for Satellites 1, 15, and 22, with 35-kHz Resolution 
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